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REMARKS 

Reconsideration of the present application is respectfully requested. Claims 
2-8, 12, 14 and 18-38 are pending. Claims 15-17, and 39 have been cancelled. 
Applicants reserve the right to pursue the content of these claims in a continuing 
application. Claims 4, 7, 8, 12, 14, 20, 24, 26, 30, 32, and 38 have been amended. 
Support for the amendments is found in the claims as originally filed, and throughout 
the specification. No new matter has been added. 

Rejections under 35 U.S.C. §112, 1^^ Paragraph - New Matter : 

Claims 24 and 32 are rejected under 35 U.S.C. §112, first paragraph as 
containing subject matter which was not described in the specification in such a way 
to reasonably convey to one skilled in the relevant art that the inventor(s) had 
possession of the claimed invention. 

Claims 24 and 32 recite a "non-human host cell". The Examiner asserts that 
"non-human" has no basis in the specification or in the claims as originally filed, and 
is considered to be new matter. 

Claims 24 and 32 have been amended to delete the term "non-human", 
therefore the rejection is obviated. 

Applicants have properly addressed by argument and amendment the 
grounds for the rejection of claims 24 and 32 under 35 U.S.C. §112, first paragraph 
and respectfully request that the rejection be withdrawn. 

Rejections under 35 U.S.C, SI 01 - Utilitv: 

Claims 2-8, 12, and 14-39 are rejected under 35 U.S.C. §101 because the 
claimed invention is not supported by either a credible or a well-established utility. 

The Examiner asserts that the claims to the isolated polynucleotide of SEQ ID 
NO: 1, a polynucleotide encoding SEQ ID NO: 2, a polynucleotide having at least 
80%, 85%, 90%, or 95% sequence identity to SEQ ID NO: 1, and a polynucleotide 
which hybridizes to SEQ ID NO: 1, wherein all the above polynucleotides encode a 
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polypeptide having ATP-dependent DNA binding activity do not nneet the utility 
requirements. Further, claims to cells, plants, and seeds comprising the above 
polynucleotides, a polynucleotide encoding at least 20 contiguous bases of SEQ ID 
NO: 2, and a polynucleotide comprising at least 30 contiguous nucleotides do not 
meet the utility requirements. The Examiner asserts that the claims do not meet the 
utility requirements for the following reasons: 

1) the predicted function is based solely upon sequence comparison with a 
RadSO gene of yeast from the prior art; 

2) no specific use of a polypeptide having ATP-dependent DNA binding activity 
or RadSO activity in a plant has been disclosed; 

3) the biological functions of a RadSO gene in plants have not been documented; 
and 

4) ATP-dependent DNA binding and RadSO activity are not specific functions. 
Applicants respectfully disagree. Applicants believe that the claims do recite 

a specific function and that this function is supported in the specification. However, 
in order to expedite examination, Applicants have cancelled claims 1S-17 and 39, 
and have amended claims 12, 14 and 20 to recite that the polypeptide is involved in 
double-strand DNA break repair. The specification at pages 1-2. especially page 1, 
lines 1S-16 and page 1, line 27 - page 2, line 2 describes the well-established utility 
of RadSO polypeptide involved in DNA double-strand break repair (see also, for 
example, abstract Ref. A2, in IDS submitted 6/23/00). RadSO is known to be 
involved in recombination as well as overall DNA repair. RadSO has also been 
shown to interact with a number of other proteins involved in DNA repair and 
homologous recombination. The present invention proposes to use the well- 
established utility of RadSO in order to modulate repair activity, and recombination, in 
order to improve gene targeting and transformation in plants (see page 2, lines 10- 
12) therefore establishing specific, substantial and credible utility for the present 
invention. 
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Applicants have disclosed a full-length Rad50 polynucleotide (SEQ ID NO: 1) 
which encodes a full-length Rad50 polypeptide (SEQ ID NO: 2) comprising known 
conserved and functional motifs (page 1-2 and Example 4). Contrary to the 
Examiner's conclusion, Applicants' identification of SEQ ID NOS: 1 and 2 as 
functional Rad50 is not based solely on one sequence comparison to one yeast 
sequence. As discussed on page 2, Applicants note that Rad50 homologues had 
been identified and characterized in mammals and in the plant Arabidopsis, As 
discussed on pages 1-2, Example 4 points out the conserved domains found in other 
RadSO proteins, and demonstrate that the identification of the SEQ ID NOS: 1 and 2 
is not merely from one sequence comparison, or one database search, but also 
involves multiple sequence comparisons and analyses, including the size of the 
encoded protein, the presence of conserved domains, the relative localization of the 
domains (/.e. positions similar to known RadSOs), percent sequence identity, and 
overall homology to known RadSO sequences. Further, in the amendment filed 
2/14/02 Applicants submitted Appendix A containing a multiple sequence alignment 
(MSA) of SEQ ID NO: 2 with other RadSO proteins from human, rat, mouse, 
Arabidopsis, yeast and C. elegans. The MSA not only showed the overall homology 
of SEQ ID NO: 2 to other known RadSO proteins, but also shows conserved amino 
acids residues and domains shared with other RadSO proteins, and also shows 
partially conserved or non-conserved amino acid residues or regions. Given the 
overall homology to other RadSO sequences and presence of conserved domains 
Applicants have disclosed and claimed sequences having credible, specific, and 
substantial utility. 

The Examiner cites a paper from Science (292:1486-87, 2001(X), copy not 
provided, previous Examiner citations from Action issued 11/19/2001 has ref (X) as 
GenBank #U668876) and cites Bork et al. (Genome Research 10:398-400 2000 (Y), 
this is ref (W) in Action issued 1 1/19/2001) to assert that the state of the art teaches 
that sequence homology alone is insufficient to detennine the functional activity of a 
gene/protein. The Examiner also asserts that further research not considered to be 
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routine would be required before one skilled in the art would know how to use the 
sequences of the claimed invention to achieve a desired trait in a plant, and that, as 
established in the courts, a utility which requires or constitutes carrying out further 
research to identify or reasonably confirm a "real world" context of use is not a 
substantial utility. 

Applicants respectfully disagree. The letter to the editor published in Science 
(Lacombe et al. 292:1486-87 2001) proposed a naming convention for plant 
glutamate receptors based on their phylogeny. Clearly, the authors do consider the 
plant sequences to be glutamate receptors. Applicants note that experimental 
evidence confirms that the plant sequences are glutamate receptor ion channels 
which conduct cations (page 1487, col. 2, and refs cited col. 3 and submitted herein 
as Appendix B). The experimental evidence functionally confirms the homology- 
based identification of the plant sequences as glutamate receptor homologues. As 
discussed in the response filed 2/14/02, Bork (Genome Res. 10:398-400, 2000) 
does not assert that homology cannot predict function, he does warn that high- 
throughput automated annotation of sequences can be inaccurate. Again, the 
Applicants stress that the identification of SEQ ID NOS: 1 and 2 is not solely based 
on an automated computer annotation, as warned by Bork. The identification of 
SEQ ID NOS: 1 and 2 as RadSO sequences is based on careful and extensive 
sequence analyses, all of which are consistent with the sequences being RadSO 
homologues. For example, SEQ ID NO: 1 encodes a full-length protein (SEQ ID 
NO: 2) with the expected molecular weight of RadSO, which further contains the 
expected sequence motifs in the expected order and location as RadSO (Walker 
boxes, nuclear localization signals, heptad repeats, and leucine zipper) and further 
has an overall sequence homology to known RadSOs as evidenced by the 
specification (e.g. pages 1-2, and Example 4) and the MSA presented in Appendix A 
submitted 2/14/02. This evidence clearly indicates that SEQ ID NOS: 1 and 2 are 
RadSO homologues containing the conserved functional domains, and meet the 
standards for utility by having credible, specific and substantial utility to use the 
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claimed RadSO sequences, having DNA repair activity, in order to modulate DNA 
repair and recombination in order to increase transformation efficiency in a plant. 
Applicants further submit that this utility extends to the breadth of the claimed 
sequences, as the claims explicitly recite the function of the encoded polypeptides 
i.e., involved in double-strand DNA break repair. Contrary to the Examiner's 
assertion, one of skill in the art will not have to conduct non-routine research in order 
to use the claimed sequences. Applicant has disclosed the sequences of SEQ ID 
NOS: 1 and 2, conserved domains and motifs (e.g., Example 4, and Appendix A), 
guidance on the construction and isolation of nucleic acids (e.g., page 10, line 29 - 
page 1 1 , line 4; page 1 3, line 30 - page 1 5, line 1 6; page 24, line 1 5 - page 26, line 
32; and page 31, line 10 - page 37, line 2), codon degeneracy (e.g., page 5, line 23 
- page 6, line 3), codon preferences (e.g., page 6, line 30 - page 7, line 1 1 ; and 
page 56, line 10 - page 57, line 2), conservatively modified variants (e.g., page 5, 
line 20 - page 6, line 25; and page 10, line 29 - page 11, line 23), conservative 
amino acid substitutions (e.g., page 6, lines 4-25; page 11, lines 6-12; page 12, line 
30 - page 13, line 2; and page 19, lines 17-33; see also Appendix A), protein 
expression (e.g., page 43, line 20 - page 47, line 25; page 49, line 28 - page 50, 
line 27; and page 53, line 1 - page 54, line 6), sequence analyses, comparisons and 
sequence identity (e.g., page 16, line 1 - page 21, line 7; page 27, lines 1-24; and 
page 57, line 28 - page 59, line 18), and vector construction, and cell and plant 
transfonnation (e.g., page 11, line 24 - page 12, line 8; page 12, lines 14-29; page 
15, lines 30-32; page 37, line 5 - page 43, line 18; page 47, line 28 - page 49, line 
26; and page 50, line 29 - page 52, line 31). As shown in an earlier response (filed 
12/17/02), assays for Rad50 were well known in the art at the time of filing (page 7- 
8). For example, assays for protein interactions (e.g. Ref. A5, A9, A10, and All, in 
IDS submitted 6/23/00) and two-hybrid screens (Ref. A8, IDS submitted 6/23/00), 
DNA binding assays (Ref. A2, IDS submitted 6/23/00), DNA repair or recombination 
assays (Ref. A1, A4, A7, and A8, IDS submitted 6/23/00), nuclease assays (Ref. 
All, IDS submitted 6/23/00), microscopy studies (Ref. A10, ID submitted 6/23/00), 
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complementation of mutants (Ref. A8, IDS submitted 6/23/00), and response to DNA 
damaging agents including chemicals and irradiation (Ref. A1, A8, and A10, IDS 
submitted 6/23/00). The guidance in the specification and the knowledge in the art 
clearly enable one of skill in the art to use the claimed sequences, relying only on 
routine screening methods. 

Applicants have properly responded by amendment and argument to the 
rejection of claims 2-8, 12, and 14-39 under 35 U.S.C. §101, and respectfully 
request that this rejection be withdrawn and not applied to pending claims 2-8, 12, 
14, and 18-38 as amended. 

Rejections under 35 U.S.C. S112, 1^* Paragraph, Enablement - Utility : 

Claims 2-8, 12, and 14-39 are rejected under 35 U.S.C. §112, first paragraph. 
Specifically, since the claimed invention is not supported by either a credible 
asserted utility, or a well established utility for the reasons set forth above, one 
skilled in the art would not know how to use the invention. 

Applicants believe they have properly addressed the rejection under 35 
U.S.C. §101 by amendment and argument and therefore respectfully request that 
the corresponding rejection of claims 12-8, 12, and 14-39 under 35 U.S.C. §112, first 
paragraph for lack of enablement be withdrawn. 

Rejections under 35 U.S.C, $112, 1^^ Paragraph. Enablement : 

Claims 2-8, 12, and 14-39 are rejected under 35 U.S.C. §112, first paragraph, 
because the specification does not enable any person skilled in the art to which it 
pertains, or with which it is most nearly connected, to make and/or use the invention 
commensurate in scope with these claims. 

The Examiner asserts that it is unclear from the disclosure how Rad50 or a 
polypeptide having ATP-dependent DNA binding activity is involved in DNA repair 
and how SEQ ID NO: 1 can be used to achieve "efficiency with which heterologous 
nucleic acids are incorporated into the genome of a target plant cell". The Examiner 
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further states that even if the Applicant shows that SEQ ID NO: 1 encoding SEQ ID 
NO: 2 has specific DNA repair activity and can be used to "modulate the efficiency 
with which heterologous nucleic acids are incorporated into the genome of a target 
plant cell" the enablement rejection will further be applied to claims broadly drawn to 
an isolated polynucleotide having at least 80%, 85%, 90% or 95% sequence identity 
to SEQ ID NO: 1, a polynucleotide that hybridizes and encodes a polypeptide having 
ATP-dependent DNA binding activity, transgenic plants, cells, seeds comprising the 
polynucleotide, a polynucleotide encoding a polypeptide comprising at least 20 
contiguous amino acids of SEQ ID NO: 2, and a polynucleotide comprising at least 
30 contiguous nucleotides of SEQ ID NO: 1 because Applicant has not provided 
guidance for how to obtain and use the polynucleotides. The Examiner asserts that 
the Applicant has not provided sufficient guidance as how to obtain any and all 
polynucleotides having the claimed structural property and still encoding a 
polypeptide having the desired functional activity. The Examiner asserts that no 
guidance has been provided for any modification to SEQ ID NO: 1 that resulted in 
the polynucleotides of claims 12, 14, 16-17, 20-21, and 39. The Examiner asserts 
that specific guidance is required for modifications to SEQ ID NO: 1 so that 
polynucleotides having both the desired structural and function characteristics can 
be obtained. 

Applicants respectfully disagree. It is clear from the specification and the art 
that Rad50 is involved in DNA repair (e.g. see pages 1-2 of the specification, and the 
IDS submitted 6/23/00, e.g. Abstract of Ref. A2). Rad50 is known to be involved in 
DNA repair and recombination, e.g. non-homologous end joining. Applicants 
propose (see page 2, lines 10-14) that by modulating the level of Rad50 protein, a 
protein involved in non-homologous end joining, one can modulate the integration of 
heterologous polynucleotides into the genome of a plant cell. Applicants have 
provided sufficient guidance for how to obtain and use the claimed RadSO 
polynucleotides and polypeptides as discussed above, including nucleic acid and 
protein isolation, determination of nucleic acid or protein sequence identity, 
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hybridizations and conserved structural features. Further, Applicants have disclosed 
the full-length polynucleotide (SEQ ID NO: 1) encoding the full-length polypeptide 
(SEQ ID NO: 2), and have deposited a plasmid comprising SEQ ID NO: 1 with the 
ATCC under Accession No. 207194. Contrary to the Examiner's assertion, one of 
skill in the art would not be required to make any and all of the claimed 
polynucleotides having the claimed structural properties and encoding a functional 
polypeptide in order to use the invention. One of skill in the art would use the 
guidance above and knowledge in the art to readily produce polynucleotides having 
at least 90% sequence identity to SEQ ID NO: 1 and encode SEQ ID NO: 2, a 
polypeptide involved in DNA repair, or readily produce polynucleotides which encode 
polypeptides having at least 90% sequence identity to SEQ ID NO: 2 and are 
involved in DNA repair. Using the disclosed sequences, the conserved domains, 
guidance in the specification, and the knowledge in the art, for example a multiple 
sequence alignment as earlier presented in Appendix A, one of skill in the art could 
identify conserved, partially conserved and non-conserved amino acid residues and 
regions, and could reasonably predict whether a conservative or non-conservative 
amino acid substitution would likely significantly impact the function of the modified 
protein compared to the unmodified protein. Only routine screening is required to 
confirm functional variants, routine screening is not undue experimentation. It is 
respectfully submitted that 35 U.S.C. §112 does not require a working example. 
Further, case law does not require that an application provide a working example. 
Nor does 35 U.S.C. §112 require that the Applicant make every embodiment. It is 
only required that the application teach the person skilled in the art how to make and 
use the invention. The present application meets that requirement. 

The Examiner asserts that the state of the prior art teaches that structural 
identity between two DNA/protein sequences does not necessarily mean that the 
sequences have the same function. The Examiner cites Lazar et al. (MCB 1988 
8(3): 1247-1 257 (U), (ref (V) in action issued 11/19/01)) and Broun et al. (Science 
1998 282:131-133 (U)) which each provide examples of very specific limited amino 
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acid changes which resulted in elimination or alteration of the experimental protein's 
catalytic activity. The Examiner also states "the nucleic acid sequences encoding 
the proteins disclosed by either Lazar or Broun would share at least 90% or 95% 
sequence identity and would hybridize to each other under the defined stringency 
conditions". Therefore, it is unpredictable whether modifications to DNA/protein will 
retain the desired functional activity and that sequence identity alone cannot be used 
to predictably determine the function of a protein/DNA. The Examiner also states 
that no transgenic plant with a desired phenotype has been disclosed. Therefore, 
given the lack of guidance, the unpredictability, lack of working examples, and the 
state of the art, one skilled would not be able to practice the invention as broadly as 
claimed. 

Applicants respectfully disagree. One of skill in the art does believe that 
structural identity, as well as the presence of functional domains and conserved 
motifs are predictive of polypeptide function, as is clearly demonstrated by the 
pervasive use of sequence searching algorithms such as BLAST, FASTA, and the 
like. As claimed, the structural identity, /.e. percent sequence identity is not the only 
criteria used, the disclosure also points out conserved functional motifs known in the 
art and the claims further recite a functional limitation that the encoded polypeptide 
is involved in DNA double-strand break repair. If the sequences of Lazar or Broun 
were similarly claimed, while the modified sequences may hybridize to the original 
sequence, they would not meet the functional limitation, and therefore would be 
explicitly excluded by the claim. Lazar discloses modifications to a human 
transforming growth factor a (TGFa) sequence, Broun discloses modifications to a 
oleate 12-desaturase sequence. Applicants wish to point out that bother references 
use the known homology to related proteins to identify and target particular amino 
acids. The references used homology to predict important conserved amino acids 
where substitution with another amino acid would likely have an impact on the 
activity of the protein. For example, Lazar's study shows that even conservative 
substitution of L48 with similar amino acids (M or I) dramatically impacted activity, as 
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predicted by the absolute conservation of leucine (L) at this position, Broun et al. 
actually note the high sequence similarity between the oleate 12-desaturase and 
oleate hydroxylase and use this to identify seven residues conserved in desaturases 
and to target them for modification of activity (see page 131 , column 2 - column 3). 
Broun et al. actually use the sequence similarity of the desaturase and the 
hydroxylase to predict which residues to change to alter the activity of the 
desaturase. In all cases, the modified protein had to be screened for the effect(s) of 
the modification. Similarly, the disclosure of SEQ ID NOS: 1-2, the conserved 
domains and motifs shown in Example 4 and known in the art, the sequence identity 
and similarity to other known RadSO sequences (for example, pages 1-2, and 
Appendix A submitted 2/14/02), the guidance on sequence analyses, comparison, 
and identity, the guidance on codon degeneracy, silent variants, and preferences, 
the guidance on conservative amino acid substitutions, the guidance on nucleic acid 
isolation and evaluation, and the ready availability in the art for assays for Rad50 
(see above) show the specification coupled with the knowledge in the art enables a 
person in the art to make and use sequences having at least 90% or 95% sequence 
identity to SEQ ID NO: 1 , or that hybridize to the full-length complement of SEQ ID 
NO: 1 and encode RadSO polypeptides involved in DNA double-strand break repair. 

The Examiner cites Amgen Inc. Chugai Pharmaceutical Co. Ltd., 18 USPQ 2d 
1016 at 1021 and 1027 (Fed. Cir. 1991). At page 1021 it is taught that a gene is not 
reduced to practice until the inventor can define it by its "physical or chemical 
properties" and at page 1027 it is taught that the disclosure of a few sequences did 
not enable claims broadly drawn to any analog thereof. 

In Amgen v. Chugai, the Federal Circuit concluded that the patent 
specification was insufTicient to enable one of ordinary skill in the art to make and 
use the invention claimed in claim 7 of the *008 patent without undue 
experimentation. As stated on page 1027, however, "it is not necessary that a 
patent applicant test all the embodiments of his invention, In re Angstadt, 537 F.2d 
498, 502, 190 USPQ 214, 218 (CCPA 1976); what is necessary is that he provide a 
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disclosure sufficient to enable one skilled in the art to carry out the invention 
commensurate with the scope of his claims. For DNA sequences, that means 
disclosing how to make and use enough sequences to justify grant of the claims 
sought." Applicants respectfully submit, that has been done in the instant 
specification. The present invention discloses how to make and use the sequences 
of the invention, /.e. sequences having at least 90% identity to SEQ ID NO: 1, as 
discussed in the paragraph above. 

The question of experimentation is a matter of degree. The fact that some 
experimentation is necessary does not preclude enablement; what is required is the 
amount of experimentation must not be unduly extensive. PPG Ina v. Guardian 
Industries Corp. (37 USPQ 1218, 1623. (Fed. Cir. 1996). The test is not merely 
quantitative, since a considerable amount of experimentation is permissible, if it is 
merely routine, or if the specification in question provides a reasonable amount of 
guidance with respect to the direction in which the experimentation should proceed 
to enable the detemiination of how to practice a desired embodiment of the invention 
claimed. Ex parte Jacl<son, 217 USPQ 804, 807 (1982 PTOBA). 

With the guidance provided in the present specification, one skilled in the art 
can readily practice the claimed invention. Therefore, it is respectfully requested 
that the rejection of claims 2-8, 12, 14-39 under 35 U.S.C. §112, first paragraph be 
withdrawn. 

Response to Arguments 

The Examiner maintains that the rejection of the claims was proper given that 
the claimed nucleic acid sequence of SEQ ID NO: 1 or nucleic acid sequences 
encoding SEQ ID NO: 2 do not have a specific utility because "ATP-dependent DNA 
binding" activity is not a specific function. Therefore, one skilled in the art cannot 
readily use SEQ ID NO: 1 or nucleic acid sequences encoding SEQ ID NO: 2 to 
achieve a desired agronomic trait in a plant without further research, therefore the 
rejection is maintained. 
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Applicants respectfully disagree with the Examiner's conclusion that "ATP- 
dependent DNA binding" is not a specific function when considered in light of the 
specification and claims directed to the RadSO sequences SEQ ID NO: 1 and 2. It is 
a specific function associated with RadSO polypeptides and is associated with the 
DNA repair activity. However, in order to expedite prosecution, Applicants have 
amended claims 12, 14 and 20 to recite that the polypeptide is involved in double- 
strand DNA break repair, a known function of RadSO. Further, Applicants believe 
they have properly addressed the rejection of the claims under 35 U.S.C 101 and 
112, 1®* paragraph, and respectfully request that these rejections be withdrawn. 

Rejections under 35 U.S.C, §112, 1^^ Paragraph, Written Description : 

Claims 2-8, 12, and 14-39 are rejected under 35 U.S.C. §112, first paragraph, 
as containing subject matter not sufficiently described in the specification to indicate 
the inventor(s) had possession of the invention. 

The Examiner asserts that the claimed invention does not meet the current 
written description requirements because ATP-dependent DNA binding activity is not 
a specific function, the specification only described SEQ ID NO: 1 the isolated 
encoding SEQ ID NO: 2, and that substantial variation is expected among 
polynucleotides comprising 30 contiguous bases of SEQ ID NO: 1 or encoding 25 
contiguous amino acids of SEQ ID NO: 2. The Examiner also states that since the 
Applicant has not described a single species of the polynucleotides of claims 12, 14, 
16-17, 20-21, and 39, the specification fails to sufficiently describe the claimed 
invention in such full, clear, concise, and exact terms that one skilled in the art would 
recognize that Applicants are in possession of the invention as broadly claimed. The 
dependent claims comprising the polynucleotides are therefore similarly not 
described. 

Applicants respectfully disagree, the specification does clearly indicate that 
the Applicants had possession of the invention as claimed. As stated above, ATP- 
dependent DNA binding activity is a specific function, however to expedite 
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examination, Applicants have amended claims 12, 14, and 20 to recite that the 
encoded polypeptides are involved in double-strand DNA break repair. Further, as 
stated above, the specification discloses full-length RadSO sequences (SEQ ID NO: 
1 & 2) comprising known functional and conserved regions, the ATCC deposited 
polynucleotide comprising SEQ ID N0:1, and further describes methods to isolate 
and analyze polynucleotides and polypeptides having a given percent sequence 
identity, and polynucleotides which hybridize to the full-length complement of SEQ 
ID NO: 1. The Examiner is reminded that every species encompassed by the 
claimed invention need not be disclosed in the specification to satisfy the written 
description requirement of 35 U.S.C. § 1 12, first paragraph. Utter v. Hiraga, 845 
F.2d 993, 6 USPQ2d 1709 (Fed. Cir. 1988). In fact, the description of a claimed 
genus can be by structure, formula, chemical name, or physical properties. See Ex 
parte Maizel, 27 USPQ2d 1662, 1669 (B.P.A.1.1992), (citing Amgen v. Chugai, 927 
F.2d 1200, 1206 (Fed. Cir. 1991)). 

Claims 12 and 17 recite that the sequence shares at least 90% sequence 
identity to the sequence of SEQ ID NO: 1 or SEQ ID NO: 2, respectively. The 
recitation of at least 90% sequence identity is a very predictable structure of the 
sequences encompassed by the claimed invention. Further, as noted above, 
Applicants have provided guidance on conserved domains, and a multiple sequence 
alignment (Appendix A, 2/14/02) of SEQ ID NO: 2 with known RadSO proteins. The 
description of a representative number of species does not require the description to 
be of such specificity that it would provide individual support for each species that 
the genus embraces. 66 Fed. Reg. 1099, 1106 (2000). Satisfactory disclosure of a 
"representative number" depends on whether one of skill in the art would recognize 
that the Applicants were in possession of the necessary common attributes or 
features of the elements possessed by the members of the genus in view of the 
species disclosed. 66 Fed. Reg. 1099, 1106(2000). Applicants submit that the 
knowledge and level of skill in the art would allow a person of ordinary skill to 
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envision the claimed invention, /.e., a sequence having at least 90% sequence 
identity to the sequence set forth in SEQ ID NO: 1 or 2. 

A genus of DMAs may be described by means of a recitation of a 
representative number of DMAs, defined by nucleotide sequence, falling within the 
scope of the genus, or by means of a recitation of structural features common to the 
genus, which features constitute a substantial portion of the genus. Regents of the 
University of California v. Eli Lilly & Co., 119 F.3d 1559, 1569 (Fed. Cir. 1997); see 
also Guidelines for Examination of Patent Applications Under the 35 U.S.C. 112, first 
paragraph, "Written Description" Requirement, 66 Fed. Reg. 1099, 1106 (2000). 
The recitation of a predictable structure of at least 90% sequence identity to SEQ ID 
NO: 1 or 2 is sufficient to satisfy the written description requirement. 

In addition, an Applicant may rely upon functional characteristics in the 
description, provided there is a correlation between the function and structure of the 
claimed invention. Id., citing Lilly at 1568. Claims 12 and 17 recite that the claimed 
sequences are, or encode, a polypeptide involved in double-strand DNA break repair 
(Rad50) thereby providing a functional characterization of the sequences claimed in 
the genus. 

Example 14 of the Revised Interim Written Description Guidelines is directed 
to a generic claim: a protein having at least 95% sequence identity to the sequence 
of SEQ ID N0:3, wherein the sequence catalyzes the reaction A — B. The Training 
Materials concludes that the generic claim of Example 14 is sufficiently described 
under § 1 12, first paragraph, because 1) "the single sequence disclosed in SEQ ID 
N0:3 is representative of the genus" and 2) the claim recites a limitation requiring 
the compound to catalyze the reaction from A B. 

Following the analysis of Example 14, Applicants submit that claims 12 and 
17 satisfy the written description requirements of § 1 12, first paragraph, one of skill 
in art would recognize that the Applicants were in possession of the necessary 
common attributes possessed by the members of the genus. The claims recite a 
defined structural parameter, i.e. percent sequence identity, and further recite a 
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functional requirement, /.e. encoding a polypeptide involved in DNA double strand 
break repair. This functional requirement is the specific description of the function of 
RadSO. Specifically, the claims of the present invention encompass sequences 
having at least 90% sequence identity to the sequence of SEQ ID NO: 1 or 2 
(structural parameter), wherein the claimed sequences encode a polypeptide 
involved double-strand DNA break repair (RadSO) activity (functional parameter). 
Claim 14 also follows the format of Example 14. The recitation of a polynucleotide 
which selectively hybridizes to the full length complement of SEQ ID NO: 1 under 
explicit high stringency hybridization and wash conditions and which encodes a 
polypeptide involved in double-strand DNA break repair also follows the guideline for 
claiming a genus of sequences by coupling a structural parameter, i.e. hybridization, 
and a functional parameter, /.e. encoding a polypeptide involved in DNA double 
strand break repair. 

Consequently, contrary to the Examiner's conclusion, the sequences 
encompassed by the genus of claims 12, 14 and 17 are defined by relevant 
identifying physical and chemical properties. In fact, the common attributes or 
features of the elements possessed by the members of the genus is that they are 
RadSO sequences and share at least 90% sequence identity to the disclosed 
sequences of SEQ ID NOS: 1 or 2, or which selectively hybridize to the full-length 
complement of SEQ ID NO: 1 under high stringency conditions. The necessary 
common features of the claimed genus are clear. 

Applicants have disclosed a full-length isolated polynucleotide and the 
encoded RadSO polypeptide in SEQ ID NOS: 1 and 2, and have provided further 
guidance on sequence isolation, analysis, and identification, codon degeneracy, 
conserved protein domains and motifs, as well as conservative amino acid 
substitutions. Applicants clearly had possession of sequences having at least 90% 
or 95% sequence identity to SEQ ID NO: 1 or 2, therefore the rejection of claims 2-8, 
12, and 14-39 are rejected under 3S U.S.C. §112, first paragraph written description 
should be withdrawn. 
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Rejections und r 35 U.S.C, $102 : 

Claim 39 is rejected under 35 U.S.C. §1 02(b) as being anticipated by 
Dolganov (in Sequence Search Result). Dolganov teaches an isolated nucleic acid 
sequence encoding a polypeptide having at least 31 contiguous amino acids of SEQ 
ID NO: 2. 

In order to expedite prosecution, claim 39 has been cancelled, thereby 
obviated the rejection under 35 U.S.C. §102. Therefore, Applicants respectfully 
request that the rejection be withdrawn. 



In light of the foregoing remarks and amendments, it is believed that claims 2- 
8, 12, 14 and 18-38 are in condition for allowance. Withdrawal of the outstanding 
rejections and allowance of all of the remaining claims is respectfully requested. 
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MolecMlar Evotatioe of Gluntamate Receptors: A Pr5miitive SigmialiEg 
Meclnamsm tlhatt Existed Before Plaints aed Animals Diverged 

Joanna Chiu,"^ Rob DeSalle,'\ Han-Ming Lam.X Lee Meisel,"^ and Gloria Coruzzi* 

♦Department of Biology, New York University; tDepartment of Entomology, American Museum of Natural History, 
New York, New York; and ^Department of Biology, The Chinese University of Hong Kong, Shatin, Hong Kong 

We performed a genealogical analysis of the ionotropic glutamate receptor (iGluR) gene family, which includes the 
animal iGluRs and the newly isolated glutamate receptor-like genes (GLR) of plants discovered in Arabidopsis. 
Distance measures firmly placed the plant GLR genes within the iGluR clade as opposed to other ion channel clades 
and indicated that iOluRs may be a primitive signaling mechanism that predated the divergence of animals and 
plants. Moreover, phylogenetic analyses using both parsimony and neighbor joining indicated that the divergence 
of animal iGluRs and plant GLR genes predated the divergence of iGluR subtypes (NMDA vs. AMPA/KA) in 
animals. By estimating the congruence of the various glutamate receptor gene regions, we showed that the different 
functional domains, including the two ligand-binding domains and the transmembrane regions, have coevolved, 
suggesting that they assembled together before plants and animals diverged. Based on residue conservation and 
divergence as well as positions of residues with respect to fimctional domains of iGluR proteins, we attempted to 
examine structure-function relationships. This analysis defined M3 as the most highly conserved transmembrane 
domain and identified potential functionally important conserved residues whose function can be examined in fiiture 
studies. 



Introduction 

Ionotropic glutamate receptors (IGIuRs) were first 
discovered in vertebrates, in which they have been 
shown to be involved in mediating fast neuronal re- 
sponses in excitatory synapses via the amino acid neu- 
rotransmitter L-glutamate (Sprengel and Seeburg 1 995). 
Two classes of animal glutamate receptors exist: iOluRs, 
which are ligand-gated ion channels, and metabotropic 
glutamate receptors (mGluRs), which are G-protein- 
linked receptors. Vertebrate iGIuRs are pharmacologi- 
cally classified into three major groups based on their 
ligand selectivity: (1) AMPA (a-amino-3-hydroxy-5- 
methyl-4-isoxazole propionate), (2) KA (kainate), and 
(3) NMDA (A^-methyl-D-aspartate). Very often, the 
AMPA and KA receptors are further grouped together 
as non-NMDA types due to their sequence similarity 
and cross-reactivity. That is, AMPA receptors can re- 
spond to kainate as well as AMPA and vice versa. Be- 
sides these three major classes of iGIuRs, there are two 
other classes of iGluRs that are related to the AMPA/ 
KA class with respect to sequence similarity: the delta 
class (Yamazaki et al. 1992; Lomeli et al. 1993) and the 
kainate-binding proteins (Gregor et al. 1989; Wada et al. 
1989). The delta class includes 81 and 62 subunits that 
were isolated in rat and mouse brain cDNA libraries. 
Although researchers have yet to demonstrate that the 
wild-type delta genes encode proteins with ion channel 
activity, a mutation in 82 shown to cause a neurodegen- 
erative disorder in Lurcher mice results in a constitu- 
tively open ion channel as tested in Xenopus laevis oo- 
cytes (Zuo et al. 1997). Using electrophysiological tech- 
niques, Rosenmund, Stem-Bach, and Stevens (1998) 
presented data that indicate iGluR channels exist as tet- 
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rameric structures, much like the voltage-gated potas- 
sium channels (Doyle et al. 1998). In contrast, other 
groups previously presented data to support a penta- 
meric structure for iGluRs (Premkumar and Auerbach 
1997) similar to the acetylcholine receptor, another li- 
gand-gated ion channel. Ionotropic glutamate receptors 
can either be homomeric or heteromeric. In both cases, 
each subunit is encoded by a single gene, and depending 
on the subunit composition of the channel, the function- 
al properties of the channel differ. All identified animal 
iGluRs contain six conserved domains that are believed 
to be functionally important. These include the two li- 
gand-binding domains (GlnHl and GlnH2, also known 
as SI and S2) that have similarity to a glutamine-bind- 
ing protein (GlnH) in Escherichia coli (Nakanishi, 
Shneider, and Axel 1990; Stem-Bach et al. 1994), and 
the four transmembrane domains (M1-M4) (fig. 1^). 
M2 is believed to be a major part of the pore and is 
important in the selective filtering of ions (Wo and Os- 
wald 1995). Unlike the other three transmembrane do- 
mains, M2 was proposed not to span the membrane, thus 
exposing the two ligand-binding domains to the extra- 
cellular side of the membrane (fig. \B) (Hollmann, Ma- 
ron, and Heinemann 1994; Bennett and Dingledine 
1995). It was suggested that the two ligand-binding do- 
mains act in concert on the extracellular side of the 
membrane to provide a ligand-binding site (Stem-Bach 
et al. 1994). 

While iGluRs were originally discovered in verte- 
brates, more recently they have been identified in in- 
vertebrates, including Caenorhabditis elegans (Maricq 
et al. 1995), and Drosophila melanogaster (Schuster et 
al. 1991; Ultsch et al. 1992, 1993). A deletion mutation 
in C. elegans glr-1 suggests a role for these receptors in 
mechanosensation and signal transduction (Maricq et al. 
1995). Surprisingly, putative iGluRs have also recently 
been identified in an organism lacking a nervous system, 
the model plant Arabidopsis thaliana (Lam et al. 1998). 
The similarity between the plant GLR genes and the an- 
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Fig. 1. — A, Diagranunatic representation of the iGluR protein showing the important functional domains, the two ligand-binding domains 
(GlnHl and GlnH2), and the four transmembrane domains (M1-M4). SP at the N-terminal end represents the signal peptide. N and C represent 
the amino- and carboxy- terminal ends, respectively. B, Diagram showing the current model of the membrane topology of animal iGluRs. C, 
Diagram of the three transmembrane domains used in the ION analysis. A Diagram showing the bacterial periplasmic amino-acid-binding 
protein and its putative region of homology with the animal iGluR proteins. 



imal iGluRs spans all the important domains that are 
conserved among the animal iGluRs, including the two 
ligand-binding domains and the four transmembrane 
segments (M1-M4). The highest similarity between 
plant GLR genes and animal iGluRs is observed in M3, 
in which the percentage of identity is above 60% (data 
not shown). Currently, there are at least four plant GLRs. 
Based on sequence similarity, the plant genes seem to 
be more closely related to the non-NMDA class of an- 
imal iGluRs, but the exact classification has yet to be 
determined functionally in a heterologous expression 
system. In planta studies using iGluR antagonists sug- 
gest plant GLRs may be involved in light signal trans- 
duction (Lam et al. 1998). The discovery of iGluRs in 
plants suggests that signaling by excitatory amino acids 
in human brains has evolved from a primitive signaling 
mechanism that existed prior to the divergence of plants 
and animals. 

In this study, we performed a phylogenetic analysis 
of the iGluR gene family of plants and animals. By in- 
cluding the plant GLR genes in the analysis, we exam- 
ined how the plant genes fit into the evolutionary history 
of these receptors. By examining the genealogy of the 
iGluR gene family, we also attempted to make predic- 
tions concerning the functional properties of the plant 
receptors based on residue conservation. Moreover, by 
looking at amino acid character state changes, we iden- 
tified residues that are invariant or are diagnostic of the 
various classes of iGluRs. Absolutely conserved resi- 
dues are likely to be functionally important, while class- 
specific variations may control subtype specificity for 
ligand. By correlating this information to the location of 
the residue within the gene as well as the functional 
properties of the various receptor classes, we identified 
possible structure-function relationships. 



Materials and Methods 

Sequences 

The animal glutamate receptor channel genes used 
in this study were obtained from Gen Bank. Table 1 lists 
all of these genes and their GenBank accession numbers. 
The four putative plant glutamate receptors were isolat- 
ed by a combination of cDNA screening and genomic 
and EST sequence analysis as described in Lam et al. 
(1998) and are included in this study. We used a two- 
step approach in examining the plant GLR genes. First, 
we aligned and compared them with amino acid se- 
quences of various kinds of ion channels (animal glu- 
tamate receptors, potassium channels, acetylcholine re- 
ceptors, and GABA^ receptors) in the database to de- 
termine if reasonable similarity to any ion channel genes 
could be established at least in the transmembrane do- 
mains (ION analysis). These genes are listed in tables 1 
and 2 with their GenBank accession numbers. The sec- 
ond step was to limit the study to glutamate receptors 
only (GLU analysis). Bacterial periplasmic amino-acid- 
binding protein gene sequences fi'om Escherichia, Sal- 
monella, and an archaebacteria (also listed in table 1) 
were obtained and used as outgroups. 

Alignment 

CLUSTAL was used to align all amino acid se- 
quences. All alignments using the default parameters in 
CLUSTAL for both the GLU and ION analyses showed 
large regions of difficult-to-align amino acids. We there- 
fore implemented a "culling'* procedure (Gatesy, De- 
Salle, and Wheeler 1993) to remove alignment-ambig- 
uous regions. 

In the ION analysis, in which we used several kinds 
of ion channel genes, only three transmembrane do- 
mains were used in our alignments in the following way. 
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Table 1 

GenBank Accession Numbers and Abbreviations for All Glutamate Receptors and Bacterial Feriplasmic Amino-Acid- 
Binding Proteins Examined in this Study 



Abbreviation Used 
Accession (if different 

Gene Species Number from gene name) 



Animal glutamate receptors 

hNMDARl Homo sapiens D 135 15 

rNMDARl Rattus norvegicus X63255 

al Mus musculus D 10028 

rNR2A /?. norvegicus M91561 

rNR2B R. norvegicus M91562 

rNRlC R. norvegicus M91563 

rNR2D R. norvegicus L31612 

dNMDARI Drosophila melanogaster X71790 

NMDAR Xenopus laevis X94156 frogglur 

rdJ R. norvegicus Z17238 

rS2 R. norvegicus Z 17239 

GluHl H. sapiens M64752 

al M. musculus X57497 

rGluRK2 R, norvegicus X54655 ratk2 

rCluRC R. norvegicus M36420 

rGluR6 R. norvegicus Z11715 

p2 M. musculus D 10054 

rKAI Rattus rattus X59996 

h-EAA2 //. sapiens S40369 humeaa2 

yl Af. musculus D 100 11 

dOluRl D. melanogaster M97192 

dGluR2 D. melanogaster M73271 

GLR 1 Caenorhabditis elegans U3466 1 celegansglr 1 

Putative GIuR C. elegans Z75545 celgansglur 

GIuR3/C Gallus gallus X89509 gallusglur3 

GluR-II Columba livia S47031 pigeonglur 

GluR4 Carassius auratus U 120 18 goldfishglur 

fGluR2Ac Oreochromis nilonica L46366 cichlidgiur 

fGluR2A Tilapia nilotica L34036 fishglur 

GluR Lymnaea stagnalis ACGAE snailglur 

GluR-Kl L stagnalis X87404 lymglur 

Plant GLR^ 

GLR I Arabidopsis thaliana AF079998 

GLR2 A. thaliana AF079999 

GLR3 A. thaliana AF007271 

GLR4 A. thaliana AC000098 

Bacterial periplasmic 
amino-acid-binding 
proteins 

glnH Escherichia coli XI 4 180 ecoliglnh 

gInP Archaeoglobus Julgidus AE001090 afggln 

glnP Salmonella typhimurium U73 1 1 1 salglnp 



The entire amino acid sequences were trimmed to in- 
clude only M1-M3 in glutamate receptors and the cor- 
responding transmembrane regions in other ion channel 
genes (fig. IQ. Since no experimental data are available 
for assigning transmembrane identity to sequences of 
the plant GLRs at the present moment, the transmem- 
brane regions of plant GLRs used in this study are pre- 
dicted based on hydropathy plots and their alignment 
with animal glutamate receptors (data not shown). These 
sequences of only approximately 1 00 amino acids were 
then aligned using CLUSTAL and a gap-to-change cost 
of 10. These alignments were of inferior quality in the 
regions of proposed transmembrane identity as deter- 
mined by CLUSTAL and showed extremely limited lev- 
els of similarity. In the GLU analysis, the entire amino 
acid sequences of all animal glutamate receptor genes 
and plant GLRs were aligned using three separate sets 
of alignment parameters (gap to change = 10, 20, 30). 



These alignments were examined for regions of align- 
ment stability (columns of amino acids that did not 
change from alignment cost to alignment cost) and am- 
biguity (as defined in Gatesy, DeSalle, and Wheeler 
1993). The alignment-ambiguous amino acid columns 
were "culled" from the data matrix. All of these prelim- 
inary alignments are available on request from the au- 
thors. 

Data Coding 

Two data-coding issues involved in the analysis of 
these sequences were the coding of outgroups and the 
coding of long gaps. A major problem we encountered 
with incorporating outgroup sequences into our data ma- 
trices was the lack of amino acid similarity in several 
regions of the genes examined. For the GLU analysis, 
we chose two bacterial and one archaebacterial peri- 
plasmic amino-acid-binding protein genes as outgroups. 
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Table 2 

GenBank Accession Numbers and Abbreviations for Non-Glutamate-Receptor Ion Channel Genes Examined in this 
Study 



Abbreviation Used 
Accession (if different 

Gene Species Number from gene name) 



Potassium channels 

MBKl (Kvl.l) Musmusculus Y00305 

AKTl Arabidopsis thaliana X62907 

ECOKCH Escherichia coli LI 2044 

dShakcr Drosophila melanogaster M 1 72 11 

KATl A. thaliana M86990 

HuK4 (Kvl.2) Homo sapiens L02752 

dShab O. melanogaster M32659 

Raw3 (rKv3.4) Rattus rattus X62841 

mShal (Kv4.1) M. musculus M64226 

MIRKl M musculus X73052 

GIRKl Rattus norvegicus L25264 

nlRKl Caenorhabditis elegans U40947 

PaKl Paramecium tetraurelia U 19907 

Acetylcholine receptors 

al Heliothis virescens AJ000399 wormal 

unc-38 C. elegans X98600 

ACR-3 C elegans Y08637 

a3 D. melanogaster Y 15593 da3 

a4-2 R, norvegicus AF007212 ra4-2 

al R, norvegicus LI 0077 ra2 

7 M musculus M30514 m7achr 

a6 Callus gallus X83889 ga6 

a8 G. gallus X52296 ga8 

S G. gallus K02903 gS 

a6 H. sapiens Y 16282 ha6 

olI H. sapiens Y08420 ha7 

GABAa receptors 

GABA-€ H. sapiens Y09765 

GABA-pl M. musculus AF024620 

GABA-p3 R, norvegicus D50671 

GABA-73 R. norvegicus X63324 

GABA-p D. melanogaster LI 7436 

GABA-rdl Aedes aegypti U28803 



Paas (1998) demonstrated that there are two major re- 
gions of supportable similarity between these bacterial 
genes and eukaryotic ionotropic glutamate receptor 
genes, and we chose to use these two regions as out- 
group sequences for rooting the eukaryotic ionotropic 
glutamate receptor analysis (fig. ID). All other regions 
of the glutamate receptor genes that were deemed align- 
ment-unambiguous (see above) that are not in these two 
regions were coded as missing in the outgroups (indi- 
cated by question marks in the matrix; fig. 2). 

In many cases, long stretches of positions in our 
alignment showed gaps. We chose to gap code these 
positions so as not to heavily weigh these regions in 
phylogenetic analysis. Our gap coding was implemented 
as unordered, and indel events were coded as distinct 
on the basis of their lengths (DeSaile and Brower 1997). 
Figure 2 shows some examples of how we implemented 
the gap coding. 

Phylogenetic Analysis 

Phylogenetic analysis using parsimony (PAUP*; 
Swofford 1998) was used to infer phylogeny. All char- 
acters were equally weighted in all analyses. This 
weighting scheme allowed us to most efficiently and 
severely test hypotheses of relationships (Kluge 1997) 



among the genes in this multiple gene family. For the 
ION analysis, uncorrected absolute distances between 
amino acid sequences of plant GLRs and various types 
of ion channel genes were calculated in PAUP* based 
on the data matrix which includes transmembrane re- 
gions only (fig. 3). Two different distance measures 
were generated in this analysis, one for M1-M3, and 
another one for M3 only. Characters 46-71 in this align- 
ment were excluded in the distance calculation from M 1 
to M3. Dglurl has an insertion in this region which cre- 
ated gaps in all other sequences. By excluding this re- 
gion, we avoid the introduction of similarity due to the 
gaps present in all of the sequences. The range and mean 
values for each comparison were then computed using 
the distances between different plant GLRs and the 
members of each ion channel class. 

For the GLU parsimony analysis, we implemented 
heuristic searches with 10 random-addition searches 
with TBR branch swapping to explore the tree space and 
narrowed the possibility that we obtained suboptimal 
trees. It was not necessary to place a limit on the number 
of trees saved in each search. When multiple equally 
parsimonious trees were obtained, we constructed strict- 
consensus trees to represent our phylogenetic hypothe- 
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Fig. 2. — Examples of gap coding. Two regions fh)m the glutamate 
receptor alignments are shown. In both examples, alignment using 
CLUSTAL resulted in the insertion of gaps longer then a single pep- 
tide. Such areas were receded so as not to outweigh these regions in 
phylogenetic analysis. A, Alignment from position 1 to position 34 
(out alignment positions) using CLUSTAL. Question marks denote 
gaps introduced by CLUSTAL. The three gaps in the plant GLR genes 
right after the space are the gaps in question. The column with 0 and 
1 shows the recoded character states for these gaps. 0 indicates that no 
gaps are present, and 1 indicates the presence of two gaps. All question 
marks arc treated as missing data in the phylogenetic analysis. B, 
Alignment from position 344 to position 368 (our alignment positions). 
Question marks denote gaps introduced by CLUSTAL. The two gaps 
in most sequences three positions to the left of the space are the gaps 
in question, and the column to the right of these two gaps represents 
the scoring system. (Accession numbers for these genes may be found 
in table 1 .) 

sis. The amino acid sequence character matrix used in 
the maximum-parsimony analysis was also analyzed us- 
ing neighbor joining (Saitou and Nei 1987). Uncorrected 
absolute distances were calculated and analyzed by 
neighbor joining in PAUP*. 

Node Robustness 

We generated three measures of node robustness 
for all nodes in our parsimony analysis: bootstrap values 
(Felsenstein 1985), Bremer (1994) indices, and jack- 
knife estimates (Farris et al. 1996). Bootstrap and jack- 
knife analyses were performed using PAUP*, and Bre- 
mer indices were calculated using AUTODECAY (Er- 
iksson 1996). Ten random-addition searches with TBR 
branch swapping were used for each replicate in boot- 
strap analysis and for the searches at each node for AU- 
TODECAY analysis. For our neighbor-joining analysis, 
robustness of nodes was inferred using 1 ,000 bootstrap 
replicates performed by PAUP*. 

Inhomogeneity Estimates 

Estimates of the congruence of different gene re- 
gions were obtained through the ILD test developed by 
Farris et al. (1994, 1995) as implemented in PAUP*. We 
partitioned our data matrix between transmembrane do- 
mains (M), GlnH domains (GlnH), and domains outside 



of these two areas (INTER) and tested for incongruence 
of these gene regions to see if there were significantly 
different phylogenetic signals emanating from the three 
different regions. 

Character Mapping 

The final data matrices were transported to 
MacClade (Maddison and Maddison 1992) and manip- 
ulated in that program. By using the "chart and trace" 
option in the program, we were able to examine the 
number of possible character states at each residue. We 
targeted regions from the transmembrane domains and 
the GlnH domains that had previously been determined 
as functional regions for analysis. In our analysis, we 
simply looked for regions of conservation as well as 
significant correlation of character changes with respect 
to changes in amino acids and functionality in the pro- 
teins. 

Results and Discussion 

Placement of Putative Plant Glutamate Receptors in 
the Context of Ion Channel Classes 

Previous distance measures demonstrated a high 
degree of identity between plant GLRs and animal i- 
GluRs, but did not accomplish a phylogenetic analysis 
(Lam et al. 1998). To determine where these plant GLR 
genes fit within the context of ion channel evolution, we 
constructed a matrix with amino acid sequences for the 
transmembrane domains from all GluRs in table 1 and 
all ion channel genes, including both ligand and voltage- 
gated ion channels, in table 2 (ION analysis). In addition 
to animal iGluRs which are shown to have extensive 
identity with the plant GLRs (Lam et al. 1998), we chose 
to examine the relationship of plant GLRs and other ion 
channel classes such as acetylcholine receptors, GABA^ 
receptors, and potassium channels due to structural sim- 
ilarities suggested in other reports (Hollmann, Maron, 
and Heinemann 1994; Bennett and Dingledine 1995; Wo 
and Oswald 1995). Glutamate receptors, nicotinic ace- 
tylcholine receptors, glycine receptors, serotonin recep- 
tors, 7-aminobutyric acid (GAB A) receptors, etc. are all 
members of the ligand-gated ion channel superfamily 
(Barnard 1992). Although in the same superfamily, i- 
GluRs are distinct in terms of primary structure (fig. 3) 
as well as membrane topology (fig. 4). The initial hy- 
pothesis based on generating hydropathy plots indicated 
that the animal glutamate receptor protein, like the nic- 
otinic acetylcholine receptors, has four transmembrane 
segments (Hollmann et al. 1989). However, more in- 
depth investigation has now led to the current model, in 
which the second transmembrane segment of iGluRs 
does not actually span the membrane (Hollmann, Ma- 
ron, and Heinemann 1994; Bennett and Dingledine 
1995). This major structural difference distinguishes the 
iGluRs from other ligand-gated ion channels, but, on the 
other hand, it suggests a closer structural and functional 
relationship between iGluRs and other ion channels. For 
example, membrane-spanning Ml and M3 together with 
non-membrane-spanning M2 in glutamate receptors fit 
the description of a structure known as a pore loop, 
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Fig. 3. — Alignment of animal and putative plant glutamate receptor transmembrane regions with several ion channel transmembrane regions 
that have putative similarity. Transmembrane regions 1-3 of glutamate, acetylcholine, and GABA^ receptors are used in this analysis. For 
potassium channels, regions S5, P, and S6 arc used due to their suggested homology with M1-M3 of animal iOluRs (Wo and Oswald 1995). 
The three transmembrane regions were aligned by determining the limits of the three regions in all genes in the figure. The three transmembrane 
regions arc labeled using the transmembrane region numbering for animal glutamate receptors. (Accession numbers for these genes may be 
found in tables 1 and 2.) 



which is common in other ion channels, such as voltage- 
gated potassium, sodium, and calcium channels, as well 
as inward rectifier potassium channels and nucleotide- 
gated cation channels (MacKinnon 1995). Wo and Os- 
wald (1995) reported conservation of amino acid resi- 
dues between glutamate receptors and potassium chan- 
nels within this functionally important region. 

The main purpose of this first analysis was to at- 
tempt to discover elements of the plant GLR transmem- 
brane sequences that might diagnose the plant GLR 
genes as glutamate receptors compared with other ion 
channels. We limited our analysis to relevant parts of 
the proteins which are possibly homologous across all 
the ion channels, i.e., the characters which make up the 
pore loop region in the potassium channels and the 
iGluRs and the corresponding transmembrane segments 
(M1-M3) in acetylcholine and GABA^ receptors (fig. 
4). The resulting alignment (fig. 3) is mistable, as evi- 
denced by drastic changes in alignment columns with 
slight changes in gap-to-change cost parameters (data 
not shown). The phylogenetic tree resulting from par- 
simony analysis is not strongly supported except for the 
node which groups the plant GLRs with the animal 
iOluRs and the nodes which group the members of the 
distinct ion channel classes (data not shown). This may 
indicate different evolutionary origins for the various 
ion channel classes, or it may simply indicate that it is 
not possible to homologize amino acid sequences in 



transmembrane domains of iOluRs and other ion chan- 
nels. In this respect, it is important to note that although 
Wo and Oswald (1995) suggested possible homology 
between transmembrane regions of potassium channels 
and animal iGluRs, our alignment leads us to suggest 
that this homology is doubtful. Instead, the similarity in 
the pore loop structure present in both classes of ion 
channels may be due to convergent evolution. 

Distance calculations of the data matrix aligned in 
figure 3 are shown in figure 4. The smaller the distance 
value, the more closely related are the sequences com- 
pared. Although the ranges for absolute distances of the 
different categories overlap slightly when considering 
M1-M3 in glutamate receptors and the corresponding 
regions in other ion channels, the mean values indicate 
that the plant GLRs are most closely related (smaller 
distance values) to animal glutamate receptors. This re- 
sult is strongly supported by previous analyses, such as 
hydropathy plots and BLAST searches (data not shown). 
This suggests that the plant GLRs have amino acid res- 
idues in their transmembrane domains that are diagnos- 
tic of animal glutamate receptor genes rather than other 
ion channel classes. In addition, when considering only 
M3 in glutamate receptors and the corresponding trans- 
membrane domain in other ion channels, the absolute 
distances between plant GLRs and animal iGluRs de- 
crease, while those between plant GLRs and non-gluta- 
mate-receptors increase dramatically. 
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Fig. 4. — Uncorrected absolute distances between plant GLRs and 
four classes of ion channels (animal iGluRs, potassium channels, ace- 
tylcholine receptors, and GABA^ receptors) based on the transmem- 
brane region alignment shown in figure 3. The membrane topology for 
each of the different ion channels in the three transmembrane segments 
considered in this analysis is also shown here. For distance calculations 
using M1-M3 of glutamate receptors and the corresponding transmem- 
brane regions in the other ion channels, characters 46-71 from the 
alignment shown in figure 3 are excluded. In this region, gaps arc 
introduced in the alignment In all sequences due to an insertion in 
dglurl (see fig. 3). For distance calculations using only M3 of gluta- 
mate receptors and the corresponding regions of other ion channels, 
characters 1-117 were excluded. In each of the comparisons, uncor- 
rected absolute distances for all different combinations of plant GLRs 
and members of each ion channel class were calculated using PAUP*. 
The mean value and the range were then computed and are presented. 
Note: The smaller the distance values, the more related the sequences. 



The alignment of the transmembrane regions of the 
various ion channels brings renewed significance to the 
similarity observed in the transmembrane regions of the 
animal iGiuRs and plant GLRs. Unlike the comparison 
between the transmembrane segments of glutamate re- 
ceptor genes and other ion channel genes, which shows 
very limited similarity, the transmembrane regions of 
animal iGluRs and plant GLRs are similar not only in 
terms of hydrophobicity, but also in terms of primary 
sequence. This evidence supports common ancestry of 
animal iGluRs and plant GLR genes rather than conver- 
gence. 

Based on the information from this analysis, we 
predict that the membrane topology of the Arabidopsis 



GLRs will be similar to that of animal iOluRs where the 
second hydrophobic segment of the protein does not 
span the membrane, thus leaving the two ligand-binding 
domains on the extracellular side of the cell membrane 
(fig. 15). This coincides with the prediction of TMPred 
analysis done on two of the plant GLR genes (Lam et 
al. 1998). Functional assays and membrane topology ex- 
periments will be required to confirm this prediction ex- 
perimentally. 

Plant Glutamate Receptor Genes in the Context of 
Glutamate Receptor Gene Evolution 

The ION analysis described above suggests that the 
plant GLR genes belong to the iGluR family. The next 
step of our analysis was to examine where they fit into 
the evolutionary history within the iGluR receptor gene 
family. We therefore constructed a matrix composed of 
amino acid sequences firom the glutamate receptor genes 
of plants and animals listed in table 1 . One of the chal- 
lenges of this analysis was to choose an appropriate out- 
group. This is not straightforward, since no putative 
iGluRs have been uncovered in the genomes of eubac- 
teria, archaebacteria, or yeast whose entire genomes 
have been completely sequenced. However, a group of 
genes known as the bacterial periplasmic amino-acid- 
binding proteins has been knovm to share identity with 
one of the functionally important modules of the animal 
iGluRs, the ligand-binding domain GlnHs (Nakanishi, 
Shneider, and Axel 1990). It was suggested that perhaps 
with the rise of multicellularity, there was increased se- 
lection pressure toward the development of membrane 
proteins necessary for cell-cell signaling (Wo and Os- 
wald 1995). We therefore used three noneukaryotic peri- 
plasmic amino-acid-binding protein sequences (two eu- 
bacteria and one archaebacteria) that share similarity 
with the putative ligand-binding domains of iGIuRs as 
outgroups (fig. ID). Visual inspection of the alignments 
we obtained using the bacterial sequences as outgroups 
reveals a much greater degree of similarity when com- 
pared with the alignment in the ION analysis (fig. 2). 
However, the regions prior to the conserved GlnHl and 
after GlnH2 showed no similarity and were easily 
"culled" (Gatesy, DeSalle, and Wheeler 1993) from the 
matrix. 

Parsimony analysis of this glutamate receptor data 
matrix using the bacterial sequences as outgroups gave 
two equally parsimonious trees of 2,686 steps with a 
consistency index of 0.705 and a retention index of 
0.782. Figure 5 shows the consensus parsimony tree 
with Bremer support and jackknife and bootstrap values. 
The plant GLR genes are consistently placed as the most 
basal group of glutamate receptors in the gene family, 
and the support for this arrangement is relatively robust 
(Bremer support at the node = 7; bootstrap at the node 
= 98%; jackknife at the node = 97%). Our neighbor- 
joining analysis generated a tree with a very similar to- 
pology and led to identical conclusion in terms of the 
placement of the plant GLRs (fig. 6). This placement 
indicates that the putative plant receptors diverged from 
the other animal iGluRs well before the various iGluRs 
that reside in animals (e.g., AMPA vs. NMDA) began 
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Fig. 5. — Phylogenetic tree generated from parsimony analysis of 35 glutamatc receptor genes and 3 bacterial periplasmic amino-acid- 
binding protein sequences. This is the consensus of two most-parsimonious trees obtained from our search. Nodes are labeled in the following 
way: Nodes with squares on them are strongly supported. The number in the square is the Bremer support value, and all nodes with squares 
had bootstrap and jackknife values of 95%- 100%. The nodes with circles on them are less strongly supported. Nodes with letters on them had 
Bremer support values, bootstrap values, and jack-knife values, respectively, as follows: A = 3/96/92, B = 2/72/77, C = 3/82/79, D = 3/85/ 
83, E = 3/67/58, F = 10/84/88, G = 1/59/58, H = 8/80/83, J = 7/86/85, K = 1/60/52, L = 9/92/96, M - 1/74/66, N = 0/50/61, P = 6/97/ 
92, and Q = 20/93/93. All other nodes with circles show the Bremer support value in the circle and at least one of the measures for robustness 
(bootstrap and jackknife) is less than 50%. 



to diverge. In addition to the placement of the plant 
GLRs, the branch lengths in the neighbor-joining tree 
(fig. 6) also indicate that whereas most of the animal 
iGluRs, with the exception of snail, C elegans, and 
Drosophila iGluRs, are very closely related, the plant 
GLRs show more significant differences. Moreover, our 
analysis establishes the relationship among members of 
the plant GLR gene family, in which GLR2 and GLR4 
are related more closely to each other than to either 
GLRl or GLRS. Members of each of the three subdivi- 
sions within the plant GLR gene family may represent 
proteins that have different functions, e.g., ligand selec- 
tivity, as in the case of animal iGluR subtypes. This 
hypothesis needs to be tested experimentally in future 
functional assays. According to this analysis, the mem- 
bers of the various classes of animal iGluRs, AMP A/ 
KA, NMDA, and delta, each reside within their appro- 
priate groups with strong node support. Previously, the 
cloned glutamatc receptor gene (glr-1) from C. elegans 
was classified into the AMPA/KA family due to se- 
quence similarity (Maricq et al. 1995). This analysis 



supports the hypothesis, although no functional data are 
available at present. 

Due to the early divergence of the plant genes, it 
is difficult to predict the function of putative glutamatc 
receptors in plants. However, based on our phylogenetic 
trees (figs. 5 and 6), it is likely that the plant GLRs will 
have ligand-binding capability, since most of the animal 
GluRs retained the ligand-binding property of the bac- 
terial periplasmic binding proteins. It is interesting to 
note that no ligand-binding or channel activity has been 
recorded with a class of animal iGluR genes known as 
the delta class. Wild-type delta genes showed no channel 
activity in heterologous systems such as X. laevis oo- 
cytes, even when glutamatc, the natural ligand for glu- 
tamatc receptors, was added to the bath solutions (Araki 
et al. 1993; Lomeli et al. 1993; Zuo et al. 1997). How- 
ever, when a mutated form of the rat 82 gene that causes 
the phenotype of the neurodegenerative Lurcher mice 
was used in similar experiments, the channels seemed 
to be constitutively open in the presence or absence of 
ligand (Zuo et al. 1997). This raises two major ques- 
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Fig. 6. — Phylogenctic tree generated from neighbor-joining analysis of 35 glutamate receptor genes and 3 bacterial periplasmic amino-acid- 
binding protein sequences showing branch lengths and bootstrap values. Uncorrected absolute distances were calculated from the same amino 
acid data matrix used in the maximum-parsimony analysis and analyzed by neighbor-joining in PAUP*. Bootstrap values arc also generated in 
PAUP* and are shown for each of the nodes except one. The node which holds all the AMPA/KA genes except for celegansglur and dglur2 
collapses in the tree from the bootstrap analysis and is not shown here. 



tions: (1) Do the plant GLRs retain ligand-binding abil- 
ity? (2) Will they respond to the same ligand as the 
animal receptors, or can they be constitutive ly open? In 
the rat 52 mutants, a single nucleotide alteration causes 
an alanine in M3 of the wild- type protein to change to 
a threonine in the mutant (Zuo et al. 1997). While most 
functional studies of animal iGluRs have concentrated 
on M2 due to the presence of RNA editing in this re- 
gion, it is significant that the greatest sequence similarity 
between animal iGluRs actually resides in M3. The 
same thing is true when the plant GLR genes are aligned 
with the animal iGluRs (fig. 4). Due to the high se- 
quence conservation of M3 among animal iOluRs and 
between plant GLRs and animal iGluRs, this region is 
suggested to have functional importance. The location 
of the mutation in the Lurcher mice (Zuo et al. 1997) 
in the M3 domain supports this hypothesis. A close ex- 
amination of the residues in the plant GLR genes cor- 
responding to the Lurcher mutation residue raises an in- 
teresting speculation. Instead of having a sequence of 
"LAA," as in wild-type delta proteins, the putative plant 
receptors have a sequence of "LTS" or "LAS" at the 



corresponding location. Like in the case of the rat 82 
mutation, which results in "LAT" at that location, the 
crucial second alanine is altered. However, instead of 
changing alanine to a threonine, as in the rat 82 mutant, 
it is substituted with serine in the plant proteins. The 
amino acid serine is very closely related to threonine in 
terms of charge and size and may confer similar func- 
tional properties. Although it is not likely that the pu- 
tative plant channels are constitutively open, as in the 
case with the rat 82 mutant, it is possible that they are 
constitutively open in certain cellular conditions, such 
as at certain extracellular ion concentrations. This and 
other possibilities are yet to be tested experimentally. 
On the other hand, perhaps these plant proteins do not 
function primarily via channel activity, but via protein- 
protein interactions. In this case, that would mean that 
channel activity actually evolved after the divergence of 
plants and animals. 

Homogeneity of the Transmembrane Domains and the 
GlnH Domains 

We tested the homogeneity of the phylogenetic sig- 
nal in the transmembrane domains versus that in the 
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Fig. 7. — A. Consensus tree generated using parsimony for the GlnH regions of 35 glutamate receptor genes. There were 16 most-parsi- 
monious trees of length 81 1; consistency index (CI) = 0.667 and retention index (RI) = 0.755. B, Consensus tree generated using parsimony 
for the regions of 35 glutamatc receptor genes that lie between the GlnH and transmembrane domains. There were eight most-parsimonious 
trees of length 1,036; CI = 0.749 and RJ = 0.810. C, Consensus tree generated using parsimony for the transmembrane regions of 35 glutamatc 
receptor genes. There were 26 most-parsimonious trees of length 699; CI = 0.735 and RI = 818. 



GlnH domains using the ILD test devised by Farris et 
al. (1994, 1995). This test was implemented to deter- 
mine if the phylogenetic signals emanating from two (or 
more) data partitions are homogeneous. We reasoned 
that an inference of homogeneity of the two kinds of 
regions would be indicative of a common evolutionary 
history or linkage. In this case, the hypothesis that these 
regions were united in an exon-shuffling scenario late in 
evolution would be rejected. In contrast, inhomogeneity 
would be indicative of incongruent evolutionary history, 
and the hypothesis that these two regions of the gene 
coevolved before plants and animals diverged would be 
rejected. Figure 7 shows the individual phylogenetic an- 
alyses of the three partitions (transmembrane domains 
[M], ligand-binding domains [GlnH], and the regions 
[INTER] between the transmembrane domains and 
GlnH domains) we erected. Visual inspection of these 
trees indicates some evident topological differences. 
However, statistical analysis of the incongruence length 
difference reveals no statistically significant incongru- 
ence among the three partitions (M vs. GlnH: P = 0.87; 
INTER vs. M: P = 0.44; INTER vs. GlnH: P = 0.84). 
These results indicate that the transmembrane domains 
and GlnH domains have most likely coevolved since 
their assembly before the divergence of plants and ani- 
mals. 

Predicting Function from Sequence 

By including the plant GLRs in our analysis, we 
have the opportunity to identify residues of the iGluR 



proteins that were conserved before the divergence of 
plant GLRs and animal iGluRs. Since these residues 
have not mutated throughout the evolution of iGluRs, 
they are most likely functionally important residues that 
are crucial to the basic operations of these receptors. The 
invariant residues in plant GLRs and animal iGluRs are 
presented in figure 8^. Conservation is apparent in a 
stretch of characters at the end of M3. Several invariant 
residues are also present in the ligand-binding domains, 
GlnHl and GlnH2. In addition, two tryptophan residues 
(W) are universally conserved in Ml and M2, respec- 
tively. Whereas at least one of these residues has been 
studied in mutagenesis experiments — mutations of R499 
result in nonfunctional channels (Uchino et al. 1992) — 
published data from mutagenesis studies are not avail- 
able on other invariant residues identified in this study. 
By conducting mutagenesis studies on these universally 
invariant residues, one should be able to gain insight 
into possible structure-fimction relationships. 

On the other hand, by recognizing the difference in 
the locations of conserved regions in the different clas- 
ses of iGluRs, we hope to identify regions of the pro- 
teins that may cause them to function differently, e.g., 
ligand and ion selectivity. In addition, by comparing the 
plant genes with the different classes of animal iGluRs 
separately (NMD A, AMPA/KA, and delta), we hope to 
gain insight into the functional properties of the putative 
plant receptors. Residues that are conserved only be- 
tween plant GLRs and each animal iGluR subclass are 
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Fig. 8. — Diagrams showing the invariant residues between plant GLRs and different animal iGluR classes: all animal iGluRs {A% NMDA 
(B), AMPA/KA (Q, and delta (D) iGluRs. In panel A, the positions of the invariant residues between plant GLRs and all animal iOtuRs are 
given, and the numbers correspond to amino acid positions in the human GluH 1 gene (GenBank accession number M64752). The approximate 
positions of the conserved residues are shown in relation to the functional domains, GlnHl, MI-M3, and GInH2, In panels B~D, residues that 
are invariant between all plant and animal genes are shown as in panel A, In addition, invariant residues between plant GLRs and each iGluR 
subclass (NMDA, AMPA/KA, and delta) are shown in bold. The residues shown in bold can be uniquely invariant between plant GLRs and 
one iGIuR subclass or invariant between plant GLRs and two of the three iGluR subclasses. By definition, these residues have to be conserved 
among members of each iGluR subclass as well. In fact, the number of residues that are conserved in each iGluR gene class is underestimated 
by adding the plant genes to each of the three analyses. The boxes representing the domains do not reflect actual sequence length. 



shown in bold in figure 85-D. The phylogenetic analy- 
ses shown in figures 5 and 6 failed to classify the plant 
GLRs into one of the gene classes of iGluRs — NMDA, 
AMPA/KA, or delta — based on a comparison using all 
of the important functional domains. However, it is still 
possible that the plant GLRs may be more closely related 
to a certain class of iGluR when looking at individual 
functional domains. According to figure 8, this is not 
the case. There appears to be more conservation be- 
tween the plant GLRs and the delta iGluRs in all func- 
tional domains. As mentioned earlier, the wild-type delta 
iOluRs have no channel activity (Araki et al. 1993; Lo- 
meli et al. 1993; Zuo et al. 1997); therefore, it is possible 
that this is the case for the plant genes. However, it is 
interesting to note that two residues, proline (P) and gly- 
cine (G), present at the end of GlnHl in NMDA and 
AMPA/KA are conserved in the plant GLRs but altered 
in the delta iGluRs. GlnHl is believed to function as 
part of the ligand-binding site in iGluRs. Since NMDA 
and AMPA/KA iGluRs are both functional in terms of 
ligand-binding, the plant GLRs may have the capacity 
for ligand binding as well, but the presence of channel 
activity that is totally identical to NMDA and AMPA/ 
KA iGluRs is another matter. It would be informative 
to conduct mutagenesis studies in which proline and gly- 
cine are altered in NMDA or AMPA/KA iGluRs to re- 
semble the delta residues, or vice versa. Such studies 
will give us answers as to whether the delta iGluRs have 
lost proper ligand-binding ability due to the loss of these 
two conserved residues, or whether residues other than 
the two mentioned are also crucial. 

Besides looking at residues in the iGluR proteins 
that are invariant between plants and animals, it is also 



interesting to look at the differences between the two. 
To examine more detailed structure-function relation- 
ships in this respect, we isolated M3 and examined it 
more closely. M3 was chosen because it is the most 
conserved region of the protein, as observed in figures 
4, 8 and 9, as well as in previous sequence alignment 
analyses (data not shown). Whereas it is believed that 
M2 serves as a major part of the channel pore (Sommer 
et al. 1991; Kohler et al. 1993), the function of M3 is 
not exactly clear yet. While most researchers have con- 
centrated their studies on M2 due to its structural im- 
portance and RNA-editing property, the analysis pre- 
sented here could identify residues that may be fonc- 
tionally important in M3. We identified residues that are 
conserved in all of the animal iGluRs (fig. 9 A), as well 
as those that are conserved in both animal iGluRs and 
plant GLRs (fig. 9B) in M3. As observed in fig, 9, the 
longest stretch of conserved residues in animal iGluRs 
is at the end of M3 (YTANLA), although the sixth res- 
idue in M3, tryptophan (W), is strikingly conserved in 
both plant GLRs and animal iGluRs. When the plant 
GLRs are added to the analysis, the conservation gen- 
erally decreases, but the level of conservation and the 
properties of the amino acids in the M3 C-terminal re- 
gion barely change. In other cases, possible amino acid 
character states added to the animal character states 
(shown in bold) by the addition of the plant genes are 
sometimes quite different in terms of amino acid charge 
and size. This is apparent in characters 1, 9, 14, 15, and 
20. For character 1 in M3, two of the four plant GLRs 
{GLRI and GLR4) have an aromatic phenylalanine (F) 
instead of a nonaromatic amino acid such as isoleucine 
(I) or valine (V). Although this character is not extreme- 
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Fig. 9. — Graphs showing the number of possible character states at each amino acid position in M3 in the GLU analysis. The number of 
possible amino acids is shown on the Y-axis, and the amino acid positions in M3 from the N-tcrminal to the C-terminal end arc shown from 
left to right on the X-axis. Character 1 in M3 corresponds to amino acid position 614 in the human GluHl gene (GenBank accession number 
M647S2). The possible amino acids at each character are shown inside the bar. The most common £imino acid at each character is located at 
the bottom of each bar, while the least common is located at the top. A, The animal graph was generated by excluding the plant and bacterial 
genes. By using the trace character function in MacClade (Maddison and Maddison 1992), each character of M3 was examined, and the number 
of possible amino acids at each position was recorded. B, The animal-and-plant graph was generated by excluding only the bacterial genes. In 
the animal-and-plant graph, the amino acids in bold arc added onto the possible character state list with the incorporation of the plant genes in 
the second graph. 



ly conserved even in animal iGIuRs with five possible 
character states, this drastic change in amino acid size 
and properties may have functional consequences. The 
case of character 9 in M3 seems to be the opposite of 
that of character 1 (fig. 9). Instead of having an aromatic 
phenylalanine (F) that is conserved in all of the animal 
iGluRs, the plant GLRs have nonaromatic leucine (L) or 
valine (V). Although these amino acids are similar in 
terms of hydrophobicity, their sizes are different. There- 
fore, if M3 is a structural part of the channel pore or 
vestibule, this character change may affect the selectiv- 
ity and permeability of ions passing through the chan- 
nels. In general, this part of the analysis points to pos- 
sible sites for mutagenesis at which residues in animal 
iGluRs are altered to resemble the corresponding plant 
residues. Results of these studies can help to uncover 
the function of residues in M3 as well as the functional 
properties of the plant GLRs. 

In conclusion, using the plant GLR genes as a tool 
to identify invariant residues may provide "fossil*' ev- 
idence with which to explore how cell-cell signaling by 
excitatory amino acids in animal brains evolved fi*om a 
primitive signaling mechanism that existed before plants 
and animals diverged. The fact that iGluRs appear to 
exist in plants and animals and not in the genomes of 
unicellular organisms whose genomes are completely 
sequenced suggests that cell-cell signaling in multicel- 



lular organisms may predate the divergence of animals 
and plants. 
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Glutamate-receptor genes in plants 



In animal brains, ionotropic glutamate 
receptors (GluRs) function as glutamate- 
activated ion channels in rapid synaptic 
transmission. We have now discovered that 
genes encoding putative ionotropic GluRs 
exist in plants, and we present preliminary 
evidence for their involvement in light- 
signal transduction. It may be that sig- 
nalling between cells by excitatory amino 
acids in animal brains evolved from a prim- 
itive signalling mechanism that existed 
before the divergence of plants and animals. 
Our findings also help to explain why neu- 
roactive compounds made by plants work 
on receptors in human brains. 

We isolated from Arabidopsis two full- 
length complementary DNAs, GLRl and 
GLR2, Each of these cDNAs encodes all of 
the signature domains of animal ionotro- 
pic GluRs, including the *three-plus-one* 



transmembrane domains (Ml to M4) and 
the putative ligand-binding domains (GbiHl 
and GlnH2), previously shown to be con- 
served between Escherichia coli glutamine 
permease (GlnH) and animal ionotropic 
GluRs* (Fig. 1). Hydropathy plots, trans- 
membrane prediction^ and protein-sorting' 
programs predict that the Arabidopsis GLR 
cDNAs encode a plasma-membrane signal 
peptide and four transmembrane domains 
(Ml to M4), of which M2 is predicted not 
to span the membrane (Fig. lb). This mem- 
brane topology is analogous to animal 
ionotropic GluRs in which the putative 
ligand-binding GlnH domains are exposed 
to the external side of the membrane^. 

In addition to their similar secondary 
structures, the Arabidopsis GLRs also share 
extensive sequence identity with animal 
ionotropic GluRs in the signature domains 



GlnHl (including the ligand-binding 
residue R), GlnH2, and Ml to M4, especially 
within M3 (Fig. la). The degree of identity 
between Arabidopsis GLRs and animal 
ionotropic GIuRs within these domains 
(63% to 16%) is similar to that between ani- 
mal ionotropic GluR subtypes that have 
kainate/AMPA rather than NMDA as ago- 
nists (60% to 32%). Arabidopsisy like ani- 
mals, seems to have a family of expressed 
genes encoding ionotropic GluRs, as judged 
by the existence of other GLR genes (Fig. la) 
and northern blot analysis (data not 
shown). GLi?- related genes also seem to 
exist in a variety of higher plants, inclu- 
ding both dicotyledons (tobacco and pea) 
and monocotyledons (rice and maize), as 
judged by heterologous Southern blot 
analysis (not shown). 

To investigate the possible in vivo 
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receptors. ^ Sequence alignment of signature 
ionotropic GluR domains conserved bet»A«en Ara- 
bidopsis GLR genes (A1GLR1-4) and animal ionotropic 
Glif? subtypes using kainate/AMPA or^JVDA as ago- 
nists. Numbers above the sequences correspond to 
residues of Arabidopsis GIBV, numbers underneath 
correspond to residues in rat GliRKI. Residues 
between the conserved domains are shown as dots. 
Arabidopsis GLRl (Genbank accession no. AF079998) 
and GLRS (AF079999) sequences are from fulHength 
cDNAs obtained by screening cDfMA libraries with 
expressed sequence tag clones (accession nos T22862 and R29880, respective- 
ly). The Gift? cDIMA corresponds to Arabidopsis genomic sequence (accession 
no. fiO0O2329). GLRS and GLR4 sequences were derived tram Arabidopsis 
genomic sequences of a bacterial and yeast artificial chromosome clones 
(AR)07271 and respectivety). Animal ionotropic Glifl genes used in 

the alignment are: AMPA/kainate subtypes, accession numbers K436418, X17184, 
X57498, M83552. X59996, S40369, U34661, M73271; or rM3A subtypes, accession 
numbers X63255. L05666. M91561. Sequences were aligned using Clustal X^. 
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h, Transmembrane organization of Arabidopsis GLRl based on hydropathy plot 
TMPred^and PSORT^ analysis. Hydrophobic domains indude a putative plasma- 
membrane signal sequence (hatched box) and four transmembrane domains 
M1-M4 (filled boxes), of which M2 is predicted not to span the membrane. Open 
boxes denote the two putative extracellular glutamate-binding domains GlnHl 
and GlnH2 (also called Si and S2), identified in animal ionotropic GluR genes as 
sharing sequence identity with Escherichia coli glutamine permease'. GInHIC 
denotes identity within the cartx»<y-terminal half of the GlnHl domain. 
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function (s) of putative ionotropic GluRs in 
plants, Arabidopsis seedlings were grown in 
the presence of DNQX, an antagonist of 
animal kainate/AMPA ionotropic GluRs. 
Plants grown on media containing DNQX 
phenocopy Arabidopsis long-hypocotyl (hy) 
mutants impaired in light-signal transduc- 
tion^, as judged by two independent crite- 
ria: DNQX at least partly blocks the ability 
of light to inhibit hypocotyl elongation (Fig. 
2a-<) and to induce chlorophyll synthesis 
(Fig. 2d-f). Plants grown in the light in the 
presence of DNQX display a significant 
light- and dose-dependent increase in 
hypocotyl elongation compared with con- 
trol untreated plants (Fig. 2b,c). This effect 
is specific for light, as the same DNQX con- 
centration has no such effect in the dark 
(Fig. 2a). The increase in hypocotyl length 
in DNQX-treated plants (1,5- to 2-fold) is 
similar to that of Arabidopsis hy mutants 



• :-43NQX' ::;vx:f^ +DNQX 




deficient in the synthesis of the phyto- 
chrome photoreceptor chromophore {hyl 
and/iy2;1.5- to 2,5-fold)''*. 

Because light and hormones interact to 
control hypocotyl elongation, one or both 
may be responsible for the DNQX effect. 
We therefore measured the effects of DNQX 
on a second light-dependent process, 
chlorophyll synthesis. DNQX treatment 
impairs the ability of light to induce the 
synthesis of chlorophyll in plants grown in 
the dark and exposed to light for 5 hours 
(Fig. 2d-f). DNQX-treated plants exhibited 
a 60% reduction in light-induced chloro- 
phyll accumulation (Fig. 2e) compared with 
controls (Fig. 2d). This effect is light spe- 
cific as growth and chlorophyll levels are 
unaffected in DNQX-treated plants grown 
in the dark (Fig, 2a,f). 

Taken together, these findings indicate 
that ionotropic GluRs may be involved in 



■ ; — DNQX (5h light) 




light-signal transduction in plants. Because 
Arabidopsis, like animals, seems to possess a 
family of GLR genes, the in planta effects of 
DNQX-treated plants may be due to inhibi- 
tion of one or more ionotropic GluR recep- 
tors in vivo. 

The existence of putative ionotropic 
GluRs in plants is surprising, and provides 
an insight into why plants make chemicals 
that act on human brains. Several classical 
ionotropic GluR agonists that activate and 
define specific ionotropic GluR subtypes 
in brains are natural plant products. 
Examples include kainate, synthesized in 
the seaweed Digenea simplex, and 
quisqualic acid, found in Quisqualis indica 
seeds^. Other less well known plant- 
derived ionotropic GluR agonists include 
p-N-methyl-amino-L-alanine, produced 
in cycads, and p-N-oxalylamino-L-alanine, 
from chick-peas, both of which are associ- 
ated with neurodegenerative diseases in 
animals^. 

It is traditionally believed that plants 
produce such neurotoxins for defence 
against herbivory. However, the discovery 
of putative ionotropic GluRs in plants and 
their possible role in light-signal transduc- 
tion suggest an alternative theory: that 
plants may synthesize ionotropic GluR ago- 
nists in order to regulate their endogenous 
ionotropic GluR receptors, and that selec- 
tive pressure in certain species then led to 
high-level production for defence against 
herbivores. 

The existence of putative glutamate 
receptors in plants raises the possibility that 
other neurologically active compounds 
made by plants (such as nicotine, cocaine 
and caffeine), which act on receptors in 
animal brains, may also act as ligands for 
endogenous receptors in plants to regulate 
a variety of as yet unknown cell-to-cell sig- 
nalling systems in higher plants, 
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Rgune 2 Effects of ionotropic Gtifl antagonist DNQX on Arabidopsis seedlings. Seedlings were sown on MS 
medium (Gibco BR_ fy/Iurashige and Skoog basal medium containing 3% sucrose and 0.9% Difco bactoagar) 
in the absence or presence of the ionotropic Glii^ antagonist CXsOX. All plates contained the so!\ent 
dimethyl sulphoxide at 16 pJ per ml medium, a-^ Effea of DNQX on hypocotyl length. Seedlings were ger- 
minated for 7 days either in the dark (a) or in a normal day- night cycle (16 h light 8 h dark) (b) at 22 in the 
presence (400 jiM) or absence of DNQX. Hypocotyl length was determined for light-grown seedlings (n, 
44-51) germinated on 0, 100. 200 and 400 jiM DNQX (c). An unpaired Rest comparing values from control 
and DNQX treatments: 100 jiM DNQX (no significant difference); 200 mM DNQX (P< 0.0001); 400 jiM DNQX 
{P< 0.0001). iM, Effect of DNQX on light-induced chlorophyll levels. Seedlings were germinated in total dark- 
ness in the absence (d) or presence (o) of 400 nM DNQX for 7 days. Plants were transferred to light for 5 h 
and chlorophyll levels quantified, in three separate pools of 30-40 seedlings grown in the dark (columns 1, 2) 
and subsequently exposed to light for 5 h (columns 3, 4) (f). An unpaired Mest indicates a significant differ- 
ence (P = 0.0031) in chlorophyll levels between DNQX-treated and untreated controls in light-treated plants, 
and no difference in dark-grown plants. 
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Eric D. Brenner, Nora Martinez-Barboza, Alexandra P. Clark, Quail S. Liang, Dennis W. Stevenson, and 
Gloria M. Coruzzi* 

Department of Biology, New York University, New York, New York 10003 (E.D.B., A.P.C., Q.S.L., G.M.C.); 
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lonotropic glutamate receptors (iCluRs) are ligand-gated ion channels that are the predominant neuroreceptors in the 
mammalian brain. Genes with high sequence similarity to animal iGluRs have been identified in Arabidopsis. To understand 
the role of Arabidopsis glutamate receptor-like (AtGLR) genes in plants, we have taken a pharmacological approach by 
examining the effects of BMAA [S(+)-)3-methyl-a, j5-diaminopropionic acid], a cycad-derived iGluR agonist, on Arabidopsis 
morphogenesis. When applied to Arabidopsis seedlings, BMAA caused a 2- to 3-fold increase in hypocotyl elongation and 
inhibited cotyledon opening during early seedling development. The effect of BMAA on hypocotyl elongation is light 
specific. Furthermore, BMAA effects on early morphogenesis of Arabidopsis can be reversed by the simultaneous applica- 
tion of glutamate, the native iGluR agonist in animals. To determine the targets of BMAA achon in Arabidopsis, a genetic 
screen was devised to isolate Arabidopsis mutants with a BMAA insensitive morphology (bim). When grown in the light on 
BMAA, bim mutants exhibited short hypocotyls compared with wild type, bim mutants were grouped into three classes 
based on their morphology when grown in the dark in the absence of BMAA. Class-I bim mutants have a normal, etiolated 
morphology, similar to wild-type plants. Class-II bim mutants have shorter hypocotyls and closed cotyledons when grown 
in the dark. Class-Ill bim mutants have short hypocotyls and open cotyledons when grown in the dark, resembling the 
previously characterized constitutively photomorphogenic mutants (cop, det, fus, and shy). Further analysis of the bim 
mutants should help define whether plant-derived iGluR agonists target glutamate receptor signaling pathways in plants. 



Glu is the predominant neurotransmitter in the 
brain. As a neurotransmitter, it activates Glu recep- 
tors at the post-synaptic membrane, which are in- 
volved in sensing environmental cues and in mem- 
ory function (Nowak et al., 1984; Isquierdo and 
Medina, 1995; Tsien et aL, 1996). Improper ionotropic 
Glu receptor (iGluR) function has been implicated in 
a variety of human diseases including Alzheimers 
and ParkinsoT\s dementia (Ikonomidou and Turski, 
1996; Forsythe and Bames-Davies, 1997). One sub- 
group of Glu receptors is comprised of the iGluRs, 
which function as Glu-gated ion channels that con- 
vey rapid synaptic transmission. iGluRs are pharma- 
cologically classified into subgroups based on agonist 
response. The two main iGluR subfamilies in animals 
are N-methyl-o-Asp (NMDA) and a-amino-3-hydroxy- 
5-methylisoxazole^-propionic acid (AMPA)/kainate 
(KA) (non-NMDA) activated iGluRs (Ikonomidou and 
Turski, 1996; Forsythe and Bames-Davies, 1997). 

In plants, it appears that Glu may also act as a 
signaling molecule. Glu supplied to plant growth 
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media has been shown to alter the expression of 
genes encoding enzymes involved in amino acid 
metabolism (Lam et al., 1994, 1998b; Oliveira and 
Coruzzi, 1999). Despite the evidence that amino acids 
may act as signals in higher plants, the mechanism of 
amino acid sensing and signaling is poorly imder- 
stood. Genes for putative amino acid sensors have 
been uncovered in Arabidopsis that have high se- 
quence similarity to ionotropic Glu receptors of ani- 
mals (Lam et al, 1998a; Chiu et al., 1999). Arabidop- 
sis GLRs have all the signature features of animal 
iGLRs, including a plasma membrane signaling pep- 
tide, two putative ligand-binding domains, and a 
"three-plus-one" transmembrane region (Lam et al., 
1998a; Chiu et al., 1999). 

To assess the function of putative Glu receptor 
genes in plants, Arabidopsis seedlings were treated 
with the iGluR antagonist 6,7 dinotropuinoxa- 
line 2,3(1H, 4H) dione (DNQX), known to block 
AMPA/KA iGLRs in animals (MuUer et al., 1989). It 
was shown that DNQX inhibits two key aspects of 
seedling photomorphogenesis in Arabidopsis: light- 
induced hypocotyl shortening and light-induced 
greening (Lam et al., 1998a). To further explore the 
targets of AtGLR function in plants, we tested 
whether other compounds known to block iGluR 
function in animals could also block aspects of Ara- 
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bidopsis growth and development. Several of these 
iGluR agonists are plant-derived products: kainate 
(KA; Monaghan et al, 1989; Bettler and MuUe, 1995) 
made by seaweed (Digenea simplex); j3-N-ox- 
alylamino-L-alanine (Ross et al., 1989) made by chick- 
peas {Lathyrus sp.), and BMAA (Copani et al., 1991) 
made by cycads. BMAA [S(+)-)3-methyl-a, )3-diamin- 
opropionic acid] has been detected in members of the 
family Cycadaceae or "cycads" (Pan et al., 1997). 
Cycads are believed to be the most primitive of gym- 
nosperms, whose remnant surviving members de- 
scended from the Mesozoic and the Paleozoic when 
cycads predominated the vegetation (Chamberlain, 
1919). BMAA was first isolated as the suspected 
cause of Parkinsonians dementia complex and amyo- 
trophic lateral sclerosis in Guam's Chamorro himian 
population, where consumption of Cycas circinalis L., 
a local food source, was prevalent (Whiting, 1963; 
Spencer et al., 1987). Subsequent to its detection and 
purification, BMAA has been shown to cause neural 
degeneration in primates when supplemented in 
their food (Spencer et al., 1987). Because BMAA is a 
natural plant product that blocks iGluR function in 
animals, we decided to test whether it would have 
any effect on plant GLRs using Arabidopsis as a 
model. In this study, we show that BMAA promotes 
hypocotyl elongation and inhibits cotyledon opening 
when applied to light-grown Arabidopsis seedlings. 
To identify the targets of BMAA action in plants, we 
used BMAA as a pharmacological tool to screen for 
Arabidopsis mutants resistant to BMAA-induced 
changes in photomorphogenesis. The isolation and 
preliminary characterization of these mutants is de- 
scribed below. 



RESULTS 

BMAA, a Cycad-Derived Glu Receptor Agonist, 
Causes a Long Hypocotyl Phenotype in Arabidopsis 

To probe the putative function of AtGLR genes in 
plants, we sought to determine whether the iGluR 
agonist, BMAA, caused any observable phenotypic 
effects on plant growth when supplied to Arabidop- 
sis seedlings. Arabidopsis seedlings were germinated 
and cultivated on Murashige and Skoog media in the 
presence or absence of BMAA. BMAA-treated seed- 
lings were evaluated for phenotypic alterations in 
treated plants compared with untreated control 
plants. At 8 d post-germination, the hypocotyls of 
seedlings grown in the light on Murashige and Skoog 
media containing 50 /xm BMAA, displayed elongated 
hypocotyls (Fig. 1 A), compared with control untreated 
plants (Fig. IB). The effect of BMAA on hypocotyl 
elongation was quantified. A dose-dependent re- 
sponse was observed at increasing concentrations of 
BMAA (Fig. 2A). In light-grown plants, a concentra- 
tion of 20 /ULM BMAA caused an increase in length of 
approximately 30%, and 50 |llm BMAA caused approx- 
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Figure 1. BMAA causes hypocotyl eiongation and inhibits cotyledon 
unfolding in light-grown Arabidopsis. Arabidopsis plants were ger- 
minated and cultivated on Murashige and Skoog media and 0.7% 
(w/v) agar containing 50 /iM BMAA (A) or no BMAA (B) for 8 d in the 
light. The average hypocotyl length for each treatment, as well as the 
arc-angle of opening for the cotyledons, is shown (n = 30). 



imately a 100% increase in hypocotyl length when 
compared with untreated plants (Fig. 2A). The effect 
on hypocotyl elongation is weaker at 100 ^iM BMAA, 
and at greater concentrations (200 /jlm) BMAA be- 
comes irJiibitory to growth (data not shown). In con- 
trast, BMAA does not induce hypocotyl elongation in 
dark-grown plants (Fig. 2B). Instead, BMAA has an 
increasingly negative effect on hypocotyl length in 
dark-grown plants at concentrations of 50 /lcm or 
greater (Fig. 2B). BMAA was also inhibitory to root 
growth (Fig. lA). BMAA also inhibits cotyledon open- 
ing in the light (Fig. lA). The arc of cotyledon opening 
is reduced to 50° in BMAA-treated plants (Fig. lA) 
when compared with 120° in control plants (Fig. IB). 

To determine whether the effects of BMAA (a Glu 
analog) could be reversed by Glu (the native agonist 
of iGluRs in animals), plants treated with 25 /lim 
BMAA were grown on media containing increasing 
amounts of l-G1u (Fig. 3). The BMAA-induced effects 
on hypocotyl elongation and inhibition of cotyledon 
opening can be partially reversed by the simulta- 
neous addition of l-G1u to the growth media (Fig. 3, 
A and B). The l-G1u reversal of the BMAA effects 
occurs in a dose-dependent manner. BMAA-induced 
hypocotyl elongation is reversed by approximately 
50% with 1 mM l-G1u and by approximately 100% 
with 10 mM L-Glu (Fig. 3A). BMAA-induced inhibi- 
tion of cotyledon opening is reversed by 20% with 1 
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Figure 2. Quantification of BMAA-induced hypocotyl elongation in 
light-grown and etiolated plants. Arabidopsis plants were grown on 
increasing concentrations of BMAA (0-100 ^im), in the light (A) or in 
continuous dark (B) for 5 d. The average hypocotyl length for each 
treatment is shown. Error bars show the se of the mean {n = 30). 



mM L-Glu, and by nearly 50% with 10 mM l-G1u (Fig. 
3B). To determine whether the reversal of BMAA 
effects on morphogenesis was specific to the amino 
acid L-Glu, we tested whether two other amino acids 
could similarly coimteract the effects of BMAA. 
L-Asp could not reverse the effects of BMAA, how- 
ever L-Gln could also reverse the effects of BMAA 
(data not shown). 

Selection of Arabidopsis Mutants Resistant to the 
Effects of BMAA on Morphogenesis 

To determine the targets and mode of action of 
BMAA in Arabidopsis, a screen was devised to iso- 
late mutants insensitive to the BMAA-induced effects 
on seedling morphology. For this screen, mutagen- 
ized (M2) Arabidopsis seeds were plated and grown 
in the light on Murashige and Skoog media contain- 
ing 50 jLLM BMAA (Fig. 4). On this BMAA-containing 
media, wild-type plants exhibit elongated hypocotyk 
and partially closed cotyledons. M2 plants that ex- 
hibited short hypocotyls and open cotyledons when 
grown on BMAA were identified as "BMAA insensi- 
tive morphology" {bim) mutants (Fig. 4A). A total of 



18,000 ethyl methanesulfonate (EMS) M2 seedlings 
were screened in the light on 50 /xm BMAA, and 10 
bim mutants were isolated (bim 18, 26, 40, 50, 59, 77, 
131, 136, 167, and 275). A representative bim mutant 
seedling, bim 26, identified in the M2 screen is shown 
in Figure 4B. 

M3 progeny from M2 bim plants were tested for 
genetic inheritance of resistance to the effects of 
BMAA. M3 progeny of two bim mutants are shown in 
Figure 5. When treated with 50 /iM BMAA and grown 
in the light, wild-type plants have elongated hypo- 
cotyls (Fig. 5A). In contrast, two representative mu- 
tants, bim 131 and bim 26, each have visibly shorter 
hypocotyls when grown in the light on 50 ijlm BMAA 
(Fig. 5A). When grown in the light minus BMAA, bim 
131 and bim 26 are indistinguishable from wild type 
as yovmg seedlings (Fig. 5B). BMAA treatment im- 
pairs root growth in wild-type plants, and in the 
majority of bim mutants. However, one mutant, bim 
50, showed partial resistance to BMAA-mediated in- 
hibition of root growth when compared with wild 
type (data not shown). 

The effect of BMAA on hypocotyl length of light- 
grown plants was quantified for all bim mutants and 
wild-type plants (Fig. 6A) and was compared with 
vmtreated plants (Fig. 6B). When grown in the light 
plus 50 ju-M BMAA; wild type has a 2- to 3-fold 
increase in hypocotyl length, compared with bim mu- 
tants (Fig. 6A). In contrast, when grown in the light 
minus BMAA, the majority of bim mutants are indis- 
tinguishable from wild type with regard to hypocotyl 
length (Fig. 6B). Only bim 131 and bim 167 have 
obviously shorter hypocotyls than wild type, when 
grown in the absence of BMAA (Fig. 6B). 



A Subset of bim Mutants Exhibit Constitutive 
Photomorphogenesis in the Dark 

The bim mutants have been grouped into three 
classes based on their dark morphology when grown 
in the absence of BMAA (Fig. 7). Hypocotyl length of 
the etiolated plants is quantified in Figure 8. Class I 
bim mutants have a normal etiolated morphology 
(elongated hypocotyls and closed cotyledons) when 
cultivated in the dark in the absence of BMAA (Fig. 
7A). This is shown for a representative bim mutant 
{bim 131) (Figs. 7A and 8A). Class-II and class-Ill bim 
mutants each have short hypocotyls when grown in 
the dark minus BMAA (Figs. 7, B and C, and 8A). 
Class-II bim mutants (bim 18, 40, 77, 136, 59, 167) have 
short hypocotyls in the dark, but their cotyledons 
remain closed (Fig. 7B). Class-Ill bim mutants {bim 26 
and 50) have short hypocotyls but also display open 
cotyledons in the dark (Figs. 7C and 8A), similar to 
the cop/det/fus mutants. The effects of BMAA on a 
representative cop mutant {copl-6), is shown in Figure 
6 and 8. BMAA has no effect on hypocotyl elongation 
of the cop 6-1 mutation in the light (Fig. 6B versus 6 A) 
or in the dark (Fig. 8B versus 8A). By contrast, BMAA 
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Figure 3. BMAA-induced hypocotyl elongation, 
as well as inhibition of cotyledon opening is 
reversed in a dose-dependent manner by GIu. 
Arabidospsis seedlings were cultivated in the 
light on Murashige and Skoog media containing 
1, 3, or 10 mM Glu in the absence {left, A and B) 
or presence (right, A and B) of 25 /llm BMAA. A, 
The average hypocotyl length for each treat- 
ment. B, The arc of cotyledon opening. Error 
bars show the se of the mean {n — 30) 
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causes slight elongation of bim mutants grown in the 
light (Fig. 6A). Thus, it appears that tihe mutation 
conferring BMAA-resistance in the bim mutants af- 
fects aspects of skotomorphogenesis. 

DISCUSSION 

We have shown that BMAA, a plant-derived agonist 
that blocks Glu receptor fxmction in animals, appears 
to alter early morphogenesis of light-grown Arabi- 
dopsis seedlings. BMAA promotes hypocotyl elonga- 
tion and inhibits cotyledon opening in the light. As 
such, BMAA-induced effects on Arabidopsis seed- 
lings phenocopy the long hypocotyl or "fty" mutants, 
defective in perceiving light and/or transmitting 
light signals in Arabidopsis (von Amim and Deng, 
1996; Fankhauser and Chory, 1997). 



We reported previously that DNQX, an antagonist 
of AMPA/KA receptors in animals, also causes a 
"/iy"-like phenotype when supplemented in the cul- 
ture media of Arabidopsis seedlings (Lam et al., 
1998a). The fact that two different iGluR interacting 
compounds (DNQX and BMAA) each induce hypo- 
cotyl elongation in light-grown seedlings provides 
support for the hypothesis that endogenous AtGLR 
genes in plants may be involved in photomorpho- 
genic development in Arabidopsis. DNQX (an iGluR 
antagonist) could potentially antagonize Arabidopsis 
GLRs, which may be involved in Hght-mediated in- 
hibition of hypocotyl growth. BMAA (an iGluR ago- 
nist) might inhibit AtGLR function but via a different 
mechanism. In animal systems, non-native agonists 
such as BMAA, can impair iGluR function because 
often iGluRs remain sensitized to these non-native 
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Figure 4. A genetic screen to isolate BMAA insensitive morphology {bim) mutants. A, Strategy to isolate Arabidopsis mutants 
resistant to the effects of BMAA is shown. Individual EMS mutagenized M2 seedlings are grown for up to 2 weeks on 
Murashige and Skoog media containing 50 fiM BMAA in the light. M2 individuals with short hypocotyls, compared with 
neighboring plants, are recovered from the BMAA-containing media and allowed to produce seed for analysis in the M3 
generation. B, Representative M2 bim plant {bim26) is shown as detected in the primary screen. 
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Figure 5. bim mutants are insensitive to BMAA-induced hypocotyl elongation. Two representative MB bim mutants (bim 13 7 
and bim 26) were cultivated for 5 d in the !ight on Murashige and Skoog media in the presence of 50 fxM BMAA (A) or in 
the absence of BMAA (B). A, Both bim 13 1 and bim 26 have a short hypocotyl phenotype (right) when compared with wild 
type (left), which has an elongated hypocotyl in the presence of BMAA. B, Plants grown in the absence of BMAA where bim 
131 and bim 26 appear similar to wild type at the early seedling stage. 



ligands (Ross et al., 1989). In contrast, animal iGluRs 
become desensitized to the native agonist, Glu. De- 
sensitization is necessary for ion channel closure and 
proper iGluR function (Geoffroy et al., 1991; Sprengel 
and Seeburg, 1995). The fact that the BMAA-induced 
effects on Arabidopsis morphogenesis are reversed 
by the addition of Glu (the native iGluR ligand), is 
consistent with the hypothesis that BMAA may act 
by blocking plant AtGLR signaling in Arabidopsis. In 
this scenario, increasing levels of Glu would compete 
with BMAA at the ligand-binding site and restore 
normal AtGLR desensitization and function. Alterna- 
tively, BMAA could act as an agonist to activate and 
open iGluR channels in plants, potentially regulating 
ion flow necessary for hypocotyl elongation. Hypo- 
cotyl expansion is largely due to increases in cell size, 
since most cells in the hypocotyl are formed during 
embryogenesis (Gendreau et aL, 1997). Previous work 
has already detected the activation of chloride chan- 
nels during hypocotyl elongation (Cho and Spalding, 
1996). In this scenario, we hypothesize that BMAA- 
induced hypocotyl elongation may be caused by ac- 
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tivation of Arabidopsis GLRs important for cell ex- 
pansion during hypocotyl elongation. 

Among other amino acids tested, we have found 
that Gin could also reverse the effects of BMAA on 
Arabidopsis growth, whereas Asp could not. Thus, 
Gin may also act similar to Glu as a potential agonist 
at a BMAA responsive site. GLn and Glu both trigger 
ion transport of Glu receptors from the cyanobacte- 
ria, Synechocystis (Chen et al., 1999), when expressed 
in a heterologous system. Glu receptors from cya- 
nobacteria show the strongest sequence similarity to 
Arabidopsis Glu receptors. Another possibility is that 
exogenously supplied Gin is assimilated and metab- 
olized to Glu, which is then able to reverse the effects 
of BMAA on Arabidopsis. In fact, HPLC analysis has 
shown that exogenously added Gin leads to signifi- 
cantly higher levels of endogenous Glu in Arabidop- 
sis (Oliveira and Coruzzi, 1999). 

To test these different hypotheses and to determine 
the targets of BMAA action in plants, we have un- 
dertaken a mutant screen in Arabidopsis using 
BMAA as a pharmacological tool. This molecular- 
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Figure 6. Quantification of hypocotyl lengths of bim mutants grown 
in the light in the presence or absence of BMAA. bim 18, 40, 77, 50, 
26, 136, 59, 131, 167, and 175 were grown for 5 d in the light in the 
presence of 50 ;i.M BMAA (A) or with no BMAA (B). Hypocotyl 
lengths of bim seedlings and wild type were measured and quanti- 
fied. Wild type is indicated with a black bar on the far right of each 
graph. cop1-6 mutant is indicated with a white bar on the far left of 
each graph. The average hypocotyl length for each treatment is 
shown (n = 30). Error bars show the se of the mean. 



during development. This phenotype is also similar 
to a number of hormone mutants deficient in growth 
including gai (Gendreau et al, 1999) and Ctrl (Kieber 
et al., 1993). The third class of bim mutants (bim 26 
and 50) has short hypocotyls and open cotyledons. 
These two characteristics are analogous to the pho- 
tomorphogenic mutants cop (Hou et al., 1993), det 
(Chory et al., 1989, 1991b), fus (Mis^ra et al., 1994; 
Kwok et al, 1996), and shy (Reed et al., 1994; Tian and 
Reed, 1999) mutants, which share these phenotypes. 
We tested the effects of BMAA on one constitutively 
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genetic approach should enable us to imderstand 
how BMAA might induce hypocotyl elongation and 
block cotyledon separation in light-grown Arabidop- 
sis. We have isolated Arabidopsis mutants insensi- 
tive to the effects of BMAA on early morphogenesis 
in the light. BMAA insensitive morphology (bfm) mu- 
tants have short hypocotyls when grown in the light 
in the presence of BMAA (Figs. 4-6). In contrast, 
wild-type plants display elongated hypocotyls under 
these conditions. The bim mutants were further sep- 
arated into three classes based on their morphology 
in the dark (Fig. 7). The first class of bim mutants has 
a normal etiolated morphology in the dark. We have 
identified two bim mutants in this class (bzm 132 and 
bim 175). The second class of bim mutant {bim 18, 40, 
59, 77, 136, 167), has short hypocotyls in the dark and 
closed cotyledons. This phenotype is similar to the 
proscute (Desnos et al., 1996) and korrigan (Nicol et al., 
1998) mutations, which affect cell wall formation 
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Dark 



Figure 7. bim mutants are subgrouped into three separate classes 
based on their dark-grown morphology, bim mutants were cultivated 
on Murashige and Skoog media in the absence of BMAA and grown 
in the dark for 5 d. Class-I bim mutants exhibit a normal etiolated 
phenotype (Fig. 7 A, right) when grown in the dark compared with 
wild type (Fig. 7 A, left). Class-I I bim mutants have a short hypocotyl 
and closed cotyledons (Fig. 7B, right) when compared with wild type 
(Fig. 7B, left). Class-Ill bim mutants have a short hypocotyl and open 
cotyledons (Fig. 7C, left) when compared with wild type (Fig. 7C, 
right). 
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Figure 8. Quantification of hypocotyl lengths of bim mutants grown 
in the dark in the absence or presence of BMAA. bim 18, 40, 77, 50, 
26, 136, 59, 131, 767, and 775 and wild type were grown for 5 d in 
the dark in the absence (Fig. 8A) or presence of 50 fXM BMAA (Fig. 
8B). Hypocotyl lengths were quantified. Wild type is indicated with 
black bars. Class-I bim mutants [bim 131 and bim 175) have wild- 
type length hypocotyls and are indicated with hatched bars. Class-ll 
and -III bim mutants, which have short hypocotyls in the dark, are 
marked with gray shaded bars. cop1-6 is shown in the far left side of 
the graph (white bar). The average hypocotyl length for each treat- 
ment is shown {n = 30), Error bars show the SE of the mean. 



photomorphogenic mutant, copl-S (Kendrick and 
Nagatani, 1991; Deng and Quail, 1992) (Figs. 6 and 
8). BMAA does not induce elongation of the copl-B 
hypocotyl. However, interpretation of these results 
must await molecular analysis of lesions in the bim 
mutants. It is important to test whether class II and 
III bim mutants are allelic to these previously char- 
acterized cop, det, fus, and shy mutants or whether 
they represent new loci. It is also important to test 
for allelism between the different bim mutants and 
to map the bim mutants to determine whether they 
are genetically linked to any AtGLR genes in 
Arabidopsis. 

An important aspect of development in seedlings 
involves the complicated interplay of light and vari- 
ous phytohormones. Auxin (Jensen et al., 1998; Kim 



et al., 1998), gibberellin (Jacobsen and Olszewski, 
1993; Steber et al, 1998), and brassinolide (Fujioka et 
al., 1997; Azpiroz et al., 1998) all act as positive 
regulators of hypocotyl elongation in Arabidopsis. 
Ethylene (Kieber et al., 1993) and cytokinins (Chory 
et al., 1991a), conversely, are believed to act as inhib- 
itors of hypocotyl elongation. Thus, the BMAA- 
mediated effects on hypocotyl length in light-grown 
plants may also involve the interaction of one or 
more of these phytohormones. It is also possible that 
BMAA blocks transduction of light signals, which 
inhibit h)rpocotyl elongation. Blue light (Liscum and 
Hangarter, 1991; Ahmad and Cashmore, 1993; Las- 
ceve, et al., 1999), red light (Somers et al., 1991; 
Nagatani et al., 1993), and far-red light (Dehesh et al., 
1993) are the major, incident wavelengths perceived 
by plants that repress hypocotyl elongation in Ara- 
bidopsis. In future studies, it will be important to 
determine whether BMAA effects are specific to one 
or more of these wavelengths in Arabidopsis. 

One interesting result from our studies is that 
BMAA induces hypocotyl elongation in Arabidopsis 
specifically in the light at low concentrations (20-50 
liM BMAA) (Fig. 2A). Because BMAA exerts its ef- 
fects at low (/LtM) concentrations, this suggests that 
BMAA could act as a signaling compound in plants. 
In species of the Cycadaceae, BMAA is detected at 
high levels (milligram BMAA /gram of tissue) (Vega 
and Bell, 1967; Duncan et al., 1989; Pan et al., 1997). 
The presence of high levels of BMAA in such tissues 
has led some researchers to suggest that BMAA may 
act as a toxin against predators (Ladd et aL, 1993). 
This theory of herbivore deterrence may explain why 
neurotoxins, such as BMAA are synthesized at high 
levels in plants. Our phylogenetic studies on GLR 
genes in plants and animals suggests that iGluRs are 
derived from a primitive signaling mechanism that 
existed before plants and animals diverged (Chen et 
al., 1999; Chiu et al., 1999). Those studies, plus the 
ones described herein, suggest the intriguing possi- 
bility that iGluR agonists made by plants may serve 
not only as herbivore deterrents, but may also act as 
signaling molecules affecting developmental pro- 
cesses in plants. We postulate, for example that 
BMAA, which appears to affect photomorphogenesis 
in Arabidopsis, may also alter light signaling in 
cycads. Whether BMAA plays a signaling role in 
cycads, or is even present at low levels in other 
species of higher plants are open questions that re- 
main to be answered. Using the Arabidopsis bim 
mutants to tmderstand the mode of action of a cycad- 
derived iGluR agonist in plants may help to address 
these questions. Furthermore, using bim mutants to 
xmderstand how BMAA mediates its effects in Ara- 
bidopsis could potentially lead to new therapeutic 
treatments of iGluR-related neurological disorders in 
htmaans. 
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MATERIALS AND METHODS 

Culture of Arabidopsis Plants 

Arabidopsis ecotype Columbia seeds were plated on 
Murashige and Skoog media (Murashige and Skoog, 
1962), 0.1% (w/v) MES [2-(N-morpholino)ethanesulfonic 
acid] (Sigma M-2933) containing 0.5% (w/v) Sue and 0.7% 
(w/v) agar. Arabidopsis seeds were placed for 2 d at 4°C 
on the growth media. Plants were grown in square (100 X 
15 mm) plates in a vertical position. Light grown plants 
were grown at 22°C during a cycle of 16-h light /8-h dark 
under cool-white fluorescent bulbs (General Electric, Fair- 
field, CT). Plants received a fluence level ranging from 40 
to 60 /lE. For dark grown seedlings, plants were incubated 
for an initial 4 h in the light to stimulate germination. 
After the light pretreatment, dark-grown plants were 
wrapped in two layers of foil and grown in the dark for 
5 d at 22°C. 

BMAA (l-BMAA hydrochloride) was purchased from 
RBI. L-Glu (Sigma G-1501), t-Asp (Sigma A-6558), and 
L-Gln (Sigma G-3126) stocks were dissolved in water and 
the pH was adjusted to 5.7. 

cop 6-1 seeds used as a control were a gift from Dr. 
Kameda, Hokkaido University (Hokkaido, Japan). 

Measurement of Hypocotyl Length and 
Cotyledon Opening 

Plants were grown in the light for 5 to 8 d. Hypocotyl 
length was measured under the view of a dissecting scope. 
The top of the hypocotyl was defined as the point where 
the petioles of the cotyledons are attached to the axis. The 
bottom of the hypocotyl was determined as the root-shoot 
jtmction (with the root being defined as the point where 
root hairs are initially apparent). Cotyledon separation was 
measured by projecting the slide image of the seedlings 
onto a screen. Two lines were then drawn along the peti- 
oles over the image of the cotyledons and the angle was 
measured where these two lines intersected (at the shoot 
apex of the seedling). 

Screen for bint Mutants in EMS Mutagenized M2 Lines 

EMS mutageruzed Columbia Co-3, glabrous seedlings 
(Lehle Seeds, Roimd Rock, TX) were cultivated as de- 
scribed above on Murashige and Skoog media (with 0.1% 
[w/v] MES, pH 5.7; 0.5% [w/v] Sue) plates containing 50 
pLM BMAA. Mutants with a BMAA insensitive morphology 
(bim) were screened after 8 to 12 d of growth, at which time 
they were transferred to Murashige and Skoog media lack- 
ing BMAA. After 1 to 2 weeks, the plants were transplanted 
to soil and allowed to set seed. 
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We have identified a homolog of the mammalian iono- 
tropic glutamate receptor genes in Arabidopsis thaliana 
{AtGluRl), This gene was found to alter Ca^^ utilization 
when overexpressed in A, thaliana. These transgenic plants 
displayed symptoms of Ca^"^ deficiency, including browning 
and death of the shoot apex, necrosis of leaf tips, and defor- 
mation of leaves. Supplementation with Ca^^ alleviated 
these phenotypes. Overall levels of Ca^^ in tissues of con- 
trol plants were not significantly different from those of 
transgenic plants, suggesting that overexpression of the 
AtGluRI gene did not affect Ca^"^ uptake. However, the rela- 
tive growth yield as a function of Ca^*^ levels revealed that 
the critical deficiency content of Ca^^ in transgenic plants 
was three times higher than that of control plants. The 
transgenic plants also exhibited hypersensitivity to Na^ and 
ionic stresses. The ion hypersensitivity was ameliorated 
by supplementation with Ca^^. The results showed that 
overexpression of the AtGluR2 gene caused reduced effi- 
ciency of Ca^"^ utilization in the transgenic plants. The pro- 
moter of the AtGluRl gene was active in vascular tissues, 
particularly in cells adjacent to the conducting vessels. This 
suggests that AtGluR2 encodes a functional channel that 
unloads Ca^^ from the xylem vessels. The results together 
suggest that appropriate expression of the AtGluR2 pro- 
tein may play critical roles in Ca^"^ nutrition by controlling 
the ion allocation among different Ca^'^ sinks both during 
normal development and during adaptation to ionic 
stresses. 

Key words: Arabidopsis thaliana — Calcium efficiency — 
Critical deficiency content — Glutamate receptor homolog — 
Ion hypersensitivity — Xylem unloading. 

Abbreviations; AtGluR2, Arabidopsis thaliana glutamate recep- 
tor homolog 2; B5 medium, Gamborg B5 medium; CDC, critical defi- 
ciency content; EST, expressed sequence tag; GUS, p-glucuronidase; 
iGIuR, ionotropic glutamate receptor, NCBI, National Center for Bio- 
technology Information. 

The nucleotide sequences reported in this paper have been sub- 
mitted to the GenBank under accession numbers AF 159498 and 
AF 159499. 



Introduction 

Calcium (Ca^^) is the second most abundant metal ion in 
plants and is an essential element for survival in higher plants. 
Millimolar concentrations of Ca^"^ are required at the outer sur- 
face of the plasma membrane to maintain its functional integ- 
rity and the structural stability of cell walls (Epstein 1972, 
Clarkson and Hanson 1980, Marschner 1995). Thus, when 
plants are subjected to Ca^^ starvation, the effects of Ca^^ defi- 
ciency are first noticed in the young and growing tissues, 
where cell division and cell enlargements are most active 
(Scaife and Turner 1984). In contrast, in the cytosol, free Ca^"^ 
is present at levels of ^10"' M or less (Kretsinger 1977). The 
cytosolic free Ca^"^ serves as an important second messenger 
that participates in signal transduction of the various environ- 
mental stresses. Numerous physiological stimuli, such as light, 
gravity, touch, cold or heat shock, oxidative stress, drought, 
osmotic shock, fungal elicitor, hormones, or salt stresses are 
known to induce transient elevations in cytoplasmic Ca^^ levels 
(Bush 1995, Sanders et al. 1999, Trewavas 1999). These tran- 
sient increases in cytosolic Ca^^ concentration are believed to 
be critical in triggering specific adaptation responses and in 
maintaining cellular homeostasis. 

Since the early 1960s, considerable research has been 
undertaken to imderstand the genetic mechanisms by which 
mineral nutrition is controlled in plants (GerlofF 1963, Epstein 
and Jefferies 1964). For example, the physiological basis 
underlying differences in Ca^"*^ nutrition has been assessed by 
screening cultivars of high Ca^"*^ efficiency from tomato, cowpea, 
and cauliflower under low Ca^*^ stress (English and Maynard 
1981, Behling et al. 1989, Horst et al. 1993, Hochmuth 1984). 
Comparisons between the inefficient and the efficient cultivars 
firom these plants provided a physiological profile associated 
with Ca^^ efficiency (Clarkson 1984, Marschner 1995). 
Although the physiological basis for Ca^"^ nutrition has been 
well described through these studies, the molecular mecha- 
nisms underlying the process of Ca^"^ uptake, root-to-shoot 
translocation, allocation of Ca^"^ among different tissues at the 
organismic level, or Ca^^ influx and its intracellular distribu- 
tion still remains elusive. 

At the organism level, Ca^"*" movement is largely depend- 
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ent on the apoplastic pathway (Kirkby and Pilbeam 1 984). Ca^^ 
is taken up mainly via the apoplastic pathway in non-differenti- 
ated root-tip tissue and during side-root emergence. The proc- 
ess by which Ca^*^ passes the Casparian band and enters the 
stele in the root zones is poorly defined. Once Ca^"^ is loaded in 
the conducting tissues, the upward movement of Ca^"^ is 
dependent on mass flow and cation exchange along the xylem 
walls. Transport of Ca^"^ is controlled by the sink strength of 
adjacent tissues, which controls the unloading of Ca^"^ firom the 
xylem sap. When the Ca^"^ supply from the roots is deficient, 
competition between Ca^"^ sinks induces Ca^^ deficiency in the 
relatively weak sinks (Clarkson 1984). Therefore it is envis- 
aged that the proper allocation among different Ca^"^ sinks is 
necessary for plant growth. Such allocation may be regulated 
by a genetically defined mechanism, although this has not yet 
been established. 

At the cellular level, Ca^"^ is passively transported into the 
cytosol through the opening of various Ca^"^ permeable chan- 
nels located at either the endo-membranes or the plasma mem- 
brane (Knight et al. 1996, Sanders et al. 1999). Families of 
Ca^^-channels have been reported through patch clamp studies, 
radioactive tracer fluxes studies, or reconstitution studies 
(White 1998, Reid and Tester 1992, Pineros and Tester 1997). 
Despite the wealth of electrophysiological information on Ca^'^- 
permeable channels, the molecular identities of Ca^'^-permea- 
ble channels in higher plants remain largely undefined. The 
wheat low-affinity cation transporter (LCTl) was shown to 
mediate Ca^^ as well as K"^ uptake in yeasts (Schachtman et al. 
1997, Clemens et al. 1998). The Arabidopsis cyclic nucleotide 
gated channels were recently identified and shown to be per- 
meable to Ca^"^ ions using human embryonic kidney (HEK) 
cells as a heterologous expression system (Kohler et al. 1999, 
Lengetal. 1999). 

There is now a sufficiently large amount of sequence data 
available for the model plant, Arabidopsis thaliana. This now 
allows the homology screening of the database as a viable 
alternative strategy for cloning plant Ca^"^ channels. Recent 
data from the Arabidopsis genome sequencing project led to 
the estimation that over 200 different plasma membrane trans- 
porters are encoded in the Arabidopsis genome (Frommer and 
Bohnert 1999). Recently, Lam et al. (1998) reported the exist- 
ence of genes encoding putative glutamate receptors in A. thal- 
iana. These genes display high sequence similarity to the iono- 
tropic glutamate receptors (iGluRs), the neurotransmitter-gated 
ion channels found in the vertebrate central nervous system. 
The iGluRs assemble to form functional channels that are 
permeable to Ca^"^ and monovalent cations when activated by 
L-glutamate (Gasic and Hollmann 1992, Wisden and Seeburg 
1993, Hollmann et al. 1994). Influx of Ca^^ through iGluRs is 
thought to play a key role in numerous physiological responses. 

In this study, we report the molecular identification and 
functional characterization of an iGluR gene homolog fi-om A. 
thaliana in the context of its physiological role in Ca^^ nutri- 
tion. 



Materials and Methods 

Genomic cloning and generation of the full-length AtGluR2 cDNA 

Arabidopsis cDNA clone 162G7T7 (partial cDNA, accession 
number R29880) was obtained fi-om the Arabidopsis Biological 
Resource Center. The identity of the clone was confirmed by sequenc- 
ing. The DNA insert (13 kb) was used to screen an Arabidopsis 
genomic library prepared in XFix. From about 1x10^ plaque-forming 
units, three independent clones were selected and purified. Bacteri- 
ophage genomic DNA was prepared by a standard purification proce- 
dure (1989). Two clones exhibited a 3.5 kb Sali fragment that 
hybridized with the cDNA probe. This fragment was cloned in pBlue- 
script KS(+) and sequenced. The full-length sequence of the cDNA 
was obtained by modified 5 '-end rapid amplification of cDNA ends 
(5'-RACE) PGR. Poly(A)^ RNA was isolated from A, thaliana as 
described (Kwak et al. 1 997a). The Marathon cDNA amplification kit 
(Clontech) was used to synthesize cDNA. The Marathon cDNA adapt- 
ers (Clontech) were ligated to the cDNA and then the resulting cDNA 
product was subjected to PCR with LA Taq polymerase (Takara) using 
the gene-specific primers, 5'-gtttca ggatccg attaagccaaagctatg-3' (primer 
I-l), and 5'-ccttggaggtccgcggaactcatcg-3' (primer 1-2), The amplified 
1 ,890 bp product was inserted into the pGEM-Teasy vector (Promega) 
and three independent clones were sequenced. The nucleotide 
sequences were compared to the genomic clone to insure that no errors 
were introduced during PCR. The fragment was digested with Bam\{\ 
and SacW, and the resulting 5'-fragment (1 to 1,853) was inserted into 
the EST clone, resulting in the full-length ^/G/m/?2 cDNA (1 to 2,975). 

Plant transformation 

Two sets of oligonucleotide primers, 5'-gtttcaggat£cgat- 
taagccaaagctatg-3' (primer 11-1), 5'-cttcccgggtctgtcatattggtctagaaggc-3' 
(primer 11-2), 5'-catccatggaacgaggccttaaaccagc-3' (primer III-l), and 
5'-taggatccccactctcccatcccatg-3' (primer III-2) were designed from the 
sequence of the AtGluR2 genomic clone. The primer set 11 was used to 
amplify 3.2 kb fragment corresponding to the AtGluR2 coding region 
from the genomic DNA. The PCR fragment was inserted into the 
pGEM-Teasy vector (Promega) and three independent clones were 
subjected to DNA sequencing analysis. This construct was digested 
with BamH\ and Sma\, and the resulting fragment was subcloned into 
the plant expression vector (pBIN-JIT), which contains a tandem 
repeat of the cauliflower mosaic virus 35 S promoter (kindly provided 
by Dr. Cathie Martin, John Innes Center, U.K.), The 1 .5 kb promoter 
region was amplified using the primer set III and the PCR product was 
digested with Nco\ and BamHl. The insert was cloned into 
pCAMBIA1303 (MRC Laboratory of Molecular Biology, U.K.). The 
resulting constructs, p35S-AtGluR2 and pl.5-GUS, were introduced 
into the AGLl strain of Agrobacterium tumefaciens and the resulting 
strains were used to transform A. thaliana ecotype Columbia (Col-0) 
by the vacuum infiltration method (Bechtold et al. 1993). Transform- 
ants were selected in Ix B5 medium (Gamborg et al. 1968) containing 
50 Jig ml"* kanamycin or 50 jig ml"' hygromicin, respectively. 

RNA extraction and RNA gel blot analysis 

Total cellular RNA was extracted from Arabidopsis seedlings 
using TRI reagent (Molecular Research Center, Inc.). Twenty micro- 
grams of RNA was separated on a denaturing 1.2% agarose gel and 
then transferred to a Hybond-N^ nylon membrane (Amersham). The 
0.3 kb fragment of the AtGluR2 gene (nucleotide 2,424 to 2,730) was 
radiolabeled with [^^PJdCTP by using a random-priming labeling kit 
(Amersham). The blot was hybridized and washed as described (Kwak 
et al. 1997b). The blot was also probed with thel8S rRNA gene of 
Brassica napus, to assess the amoimt of RNA. 



76 AtGluR2, an Arabidopsis glutamate receptor gene 



Table 1 Putative glutamate receptors in Arabidopsis genome 
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Plant growth 

Plants were either grown in a composite soil in a greenhouse, or 
grown on sterile agar plates in petri dishes (22°C/19°C, 16 h light/8 h 
dark cycle, 100 jimol m~^ s~' photon flux density). For studies of Ca^^ 
supplementation, a soil-saturating volume of 4 mM CaClj was added 
every 2 d. For ionic sensitivity analysis, surface-sterilized seeds were 
grown for 2 weeks in Ca^^-depleted l/4x B5 medium, supplemented 
with the appropriate ions. 

Growth measurement and C<i^^ analysis 

Seedlings were grown on sterile agar plates as described above. 
The aerial portions of the plants were harvested and the fresh weight 
was measured individually. The tissues were collected and dried at 
70°C for 72 h. The dry weight was determined and then the dried tis- 
sues were digested with 65% nitric acid at 200X. The resuhing diges- 
tate was resuspended in 1 M nitric acid. The total Ca^* content per 
gram of dry weight was determined by an atomic absorption spectro- 
photometer (Unicam Atomic Absorption, U.K.). 

GVS assays 

GUS histochemical staining was performed on whole seedlings. 
Tissues were vacuum infiltrated in a buflFer containing 1 mg ml'* X- 



Gluc (5-bromo-4-chloro-3-indolyl-p-gIucuronide), lOOmM sodium 
phosphate pH7, 1% Triton XI 00, 5 mM potassium ferrocyanide, 
5 mM potassium ferricyanide, and 1 mM EDTA. The enzymatic reac- 
tion was performed overnight at 37**C in the dark. After washing in a 
buffer containing 100 mM sodium phosphate (pH 7) and 1 mM EDTA, 
tissues were cleared through a graded ethanol series. Sections (7 ^m 
thick) of GUS-stained material were prepared on paraffin embedded 
tissues and cut with an LKB microtome. 

Results 

Families of ionotropic glutamate receptor homologs from A. 
thaliana 

In a previous study from our laboratory, a plant homolog 
of mammalian ionotropic glutamate receptors was isolated 
from Brassica napus (Kwak et al. 1997a). The Brassica 
expressed sequence tag (EST) was used as a query to search for 
related sequences from the GenBank database of A. thaliana. 
This search unveiled one EST (accession no R29880) and 19 
genomic sequences for iGIuR homologs (Table 1). 
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Molecular cloning of an Arabidopsis homolog of ionotropic 
glutamate receptors 

The identified Arabidopsis EST corresponded to a 1 .3 kb 
fragment of the C-terminal sequence. This DNA fragment was 
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used to screen both a genomic library and a cDNA library. A 
genomic clone (accession no AF 159499) was obtained, but no 
hybridization was obtained in screening of the cDNA library. 
Thus, the full-length cDNA sequence was obtained by 5'- 
RACE-PCR. The full-length cDNA (accession no AF 159498, 
nucleotide 1 to 2,975) was reconstituted from the cDNA and 
PCR clone. Sequence comparison of the genomic clone and the 
cDNA was performed to verify the accuracy of the PCR prod- 
uct. The sequence contained a deduced open reading frame 
encoding a polypeptide of 912 amino acid residues with a pre- 
dicted molecular mass 102 kDa (Fig. lA). Comparison of the 
predicted amino acid sequence of AtGluR2 against the NCBI 
non-redundant protein database revealed that the AtGluR2 
sequence shares low (21% identity) but significant similarity to 
that of manmialian iGluRs throughout the transmembrane 
regions (62% identity in the TM III region) and pore-forming 
region of the peptide (Fig. IB). The genomic organization of 
the AtGluR2 suggested that this gene also displays a distinc- 
tive exon-intron structure similar to that of the kainate sub- 
types in the pore-forming domain (Fig. IC; Sommer et al. 
1991, Herb et al. 1996). When hydropathy plots were con- 
structed from the deduced amino acid sequence, AtGluR2 was 
predicted to display a three-plus-one transmembrane structure, 
as shown in Fig. ID. This structure is consistent with the struc- 
ture expected for ionotropic glutamate receptors (Hollmann et 
al. 1994, MacKinnon 1995). 



Fig. 1 Deduced amino acid sequence and the hydropathy analysis of 
AtGliiR2, and the genomic structure of the AtGluRl gene. (A) The 
deduced amino acid sequence of AtGluR2 (GenBank accession no 
AF 159498). The sequences corresponding to putative transmembrane 
domains (TM I, TM III, and TM IV) are underlined. The putative 
pore-forming domain (Pore) at residues 620-640 is double-under- 
lined. The arrowheads indicate the positions of the exon-intron junc- 
tions identified in the genomic clone by comparison with the cDNA 
sequence. (B) Sequence alignment of the pore-forming region and the 
transmembrane domain III among plant glutamate receptor homologs 
and mammalian iOluRs. Numbers above the sequence correspond to 
amino acid residues of AtGluR2. Residues between the two domains 
are shown as dots. Identical and similar residues are indicated by dark 
and light shades, respectively. Arabidopsis glutamate receptor 
homologs used in the alignment are AtGluR2, GLR2, GLR6, and 
GLUR3 (GenBank accession nos AF159498, AF079999, AF170494, 
and AFl 67355). BnGluRl and NtGluRl sequences were derived from 
a Brassica napus cDNA (accession no AFl 09392) and a tobacco par- 
tial cDNA (unpublished data) clone. Mammalian iGluRs are AMPA/ 
kainate subtypes, GluRl, GluR2, and GluR5 (GenBank accession nos 
X17184, M85035, and M83560). (C) Schematic diagram of the 
genomic clone (GenBank accession no AFl 59499). The boxs and the 
bold lines indicate the exons and the introns, respectively. The puta- 
tive transmembrane domains are marked with dots. The pore-forming 
domain is underlined. (D) The hydropathy analysis of the AtGluR2 
peptide sequence. The hydropathy plot is generated with TMpred 
(http://www.ch.embnet.org^sofhvare/TMPRED_form.html). The dots 
and the arrowhead indicate the putative transmembrane domains and 
the pore-forming domain, respectively. 
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Fig. 2 Expression of AtGluR2 in transgenic Arabidopsis plants. Total 
RNA was extracted from 3-week-old seedlings of seven representa- 
tive lines and 20 \i% of total RNA were analyzed for the AtGluR2 gene 
expression. The RNA gel blot was hybridized with 0.4 kb fragment of 
the AtGluR2 gene (nucleotide 2,424 to 2,730). Endogenous AtGluR2 
transcripts were not detected in wild-type Colimibia plants and trans- 
genic lines expressing the vector (control) alone. 



AtGluR2 expression levels in transgenic Arabidopsis 

To evaluate the physiological roles of the AtGluR2 gene in 
plants, we generated transgenic Arabidopsis lines that stably 
express the sense-oriented AtGluR2 cDNA under the control of 
the dual 35S promoter of the cauliflower mosaic virus. During 
the initial screening on selective media, 1 9 independent trans- 
genic lines displayed distinctively abnormal development. 
Seven representative transgenic lines were tested for the 
AtGluR2 expression by performing RNA gel blot analysis (Fig. 
2, upper row). In wild-type Arabidopsis and the vector-trans- 
formed control plants, the AtGluR2 mRNA expression could 
not be detected from total RNA extracted from whole seed- 
lings. The homozygous transgenic lines showed increased 
expression of the AtGluR2 transcript, although the level of 
expressions varied across the transgenic lines. Thus, the 
AtGluR2 gene was successfully overexpresed in trangenic Ara- 
bidopsis lines. 

AtGluR2 overexpression disrupts normal development 

Figure 3 shows the growth characteristics of two repre- 
sentative transgenic lines. Transgenic plants exhibited symp- 
toms that are highly consistent with Ca^^ deficiency. When 
plants are subjected to Ca^"^ deficiency, the symptoms are easily 
noticed in the growing points of the upper parts of plants. 
Brown melanin compounds resulting from polyphenolic 
oxidation are also found in the Ca^"^ deficient tissues (Kirkby 
and Pilbeam 1984). As shown in Fig. 3B, the G-5 line dis- 
played browning and necrosis of the growing points, curling 
and deformation of leaves, and disintegration of petioles, when 
seedlings were grown on Ix 35 medium for 2 weeks. When 
this line was grown in composite soil for 4 weeks, it produced 
necrotic and deformed leaves and many lateral buds. The 
whole growth was also severely retarded (Fig. 3C). The A- 12 



line developed milder symptoms than the G-5 lines, exhibiting 
bushy and dwarf stature with large numbers of short secondary 
inflorescences. Necrosis at the tips and margins of young 
leaves was also apparent (Fig. 3D). The fertility of transgenic 
lines was also reduced and this was positively correlated with 
the severity of other morphological defects (data not shown). 

Ca?^ supplementation mitigates AtGluR2-induced symptoms 

Since the morphological phenotypes of the transgenic 
lines was similar to the Ca^^ deficiency symptoms, we tested 
whether providing Ca^^ as a supplement could alleviate the 
morphological alterations caused by AtGluR2 overexpression 
in the transgenic plants. After seed germination, plants were 
watered in the presence or absence of 4 mM CaClj every 2 d. 
This concentration of Ca^^ is equivalent to that in Hoagland*s 
nutrient solution (Hoagland and Amon 1950). Transgenic A- 12 
and C-1 plants supplemented with Ca^^ displayed developmen- 
tal patterns similar to that of plants carrying the empty vector, 
although the extent of amelioration provided by Ca^"^ supple- 
mentation varied across the individual lines (Fig. 4). Addition 
of Ca^"^ mitigated all abnormalities seen in the transgenic 
plants, including those in necrosis and deformity of leaves, 
repressed growth of stem elongation, reduced apical dominance 
and lowered fertility. The plants carrying the empty vector 
showed similar growth characteristics in the presence or 
absence of Ca^"^ supplementation, although with extra supply of 
Ca^"*" the plants grew more vigorously. The result confirmed 
that the symptoms exhibited by overexpression of the AtGluR2 
gene are due to Ca^"^ deficiency. 

TTie Cc^^ levels in shoots were similar in control and transgenic 
plants 

To test whether the Ca^^ deficiency seen in the transgenic 
lines are due to disturbance Ca^"^ accumulation, we measured 
Ca^^ contents as a function of Ca^^ concentration in the grovrth 
media. Plants were grown for two weeks in Ca^^ depleted l/4x 
B5 medium and then the medium was supplemented with 
10 ^M, 100 ^iM, 250 ^M, 1 mM and 6 mM of CaClj, respec- 
tively. The control plants did not develop any signs of Ca^^ 
starvation at any of the Ca^^ concentrations tested. In contrast, 
the transgenic plants displayed Ca^^ deficiency symptoms 
when the external Ca^*^ concentration was below 100 ^iM. 
Since the symptoms of Ca^"^ deficiency were primarily 
observed in aerial parts of the plants, the chemical analysis was 
performed with samples from the shoots. As shown in Fig. 5A, 
control and transgenic plants had little differences in Ca^^ lev- 
els in the shoots. Therefore, the differential responses to low 
Ca^"^ nutrition between the control and transgenic plants could 
not be explained by differences in Ca^"^ within the plant shoots. 
This indicates that the Ca^^ deficiency phenotyj>e in the shoot 
of the transgenic plants was not attributable to reduced accu- 
mulation or uptake of Ca^^. 
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Fig. 3 Phenotypes of transgenic Arabidopsis overexpressing the AtGluR2 gene. (A) The control plant transformed with the vector alone. (B) 
Seedling of the transgenic line G-5 grown on 1 x BS medium for 2 weeks. Necrosis and browning of the shoot apex and collapse of petioles were 
observed. Bar: 1 mm, (C) Close-up view of the soil-grown transgenic line G-5. Curled and deformed leaves and many undeveloped lateral buds 
were present. Bar: 1 cm. (D) The transgenic line A- 12. Necrosis was evident at the tips and margins of young leaves. The plant exhibited a stunted 
and bushy stature with large numbers of short secondary inflorescences. 
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Fig. 4 Supplementary Ca^^ ameliorates the AtGluRl-iMuctd symp- 
toms. The vector control and the transgenic A- 12 and C-1 plants over- 
expressing the AtGluR2 gene were grown in a composite soil in a 
greenhouse for 3 weeks. Plants were watered every 2 d in the presence 
or absence of 4 mM Ca^^. 



The critical deficiency content (CDC) is higher in transgenic 
plants than that of control plants 

To further reveal the physiological basis of the Ca^^ defi- 
ciency symptoms seen in the transgenic lines, we measured the 
degree of the efficiency of Ca^^ utilization in the plant shoots 
(Gerlofif and Gabelman 1983, Marschner 1995). For this pur- 
pose, we characterized the relationship between Ca^"^ content in 
shoots and relative shoot growth. When 90% of the maximal 
growth was taken as a reference point, the CDC of Ca^^ for the 
control plants and transgenic A- 12 and C-1 plants were approx- 
imately 4 mg g'' dry weight, 13 mg g~* dry weight, and 12 mg 
g"' dry weight, respectively (Fig. 5B). Thus, the CDC is three- 
fold higher in the transgenic plants. This result together with 
the data in Fig. 5 A demonstrates that overexpression of the 
AtGIuR2 gene induces Ca^"^ deficiency symptoms by impair- 
ing the efficiency of Ca^"^ utilization in the transgenic plant. 

Overexpression of AtGluRl renders transgenic plants 
hypersensitive to IC and Na^ salt 

Ca^^ starvation in plants is rarely caused by a complete 
absence of Ca^"^ in the soil. Rather, the deficiency is more often 
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Fig. 5 The effect of Ca^^ concentration on plant growth. Plants were 
grown in the l/4x B5 medium supplemented with various concentra- 
tions of CaClj. After 2 weeks of growth, the fresh weight was meas- 
ured individually to analyze the relative growth. Seedlings were then 
pooled to yield the required dry weights for the chemical analysis. 
Ca^^ contents were determined with an atomic absorption spectropho- 
tometer, (A) The effect of external Ca^^ on the concentration of Ca^^ in 
the aerial parts. Shoots were collected by excision below the rosette 
from whole seedling plants. Data represent the averagelSD of three 
independent assays. (B) The relative growth as a function of Ca^* con- 
centration in the aerial parts. Data are expressed as the mean±SE of 
three independent experiments (25 seedlings were used for each meas- 
urement). Shown are values relative to the fresh weights of seedlings 
grown on l/4x B5 medium supplemented with 6 mM CaClj. 



induced by low concentrations of Ca^^ in relation to other cati- 
ons, such as high K"^, Na^, or Mg^^ (Adams and Ho 1990, Mass 
and Grieve 1987, Alcantara and De La Guardia 1987, Hirschi 
1999). Increase of the external concentration of these elements 
of^en reduce the absorption of Ca^"^ and thereby induce Ca^* 
deficiency in salt-stressed plants (Robson and Pitman 1983). 
Furthermore, Ca^"^ deficiency is often diagnosed by an 
increased sensitivity to other cations in crop plants (Scaife and 
Turner 1984). We, thus, tested sensitivity of the transgenic lines 
to other cations. The transgenic plants were grown in l/4x B5 



medium with addition of K"^, Na"^, and Mg^"^ for 2 weeks. The 
fresh weight of seedlings grown in l/4x B5 medium was simi- 
lar in the control and transgenic plants. The relative growth 
was, then, analyzed as a function of concentration of these cati- 
ons. The results showed that the transgenic plants showed nota- 
ble hypersensitivity to and Na^ (Fig. 6A, B). The sensitivity 
of the transgenic plants to Mg^*^ was not significantly distin- 
guishable from the control plants (Fig. 6C). The relatively dif- 
ferent interaction of the Ca^^ with other cations has been 
observed in other plants, too (Hirschi 1999). The AtGluR2 
transgenic plants were not notably hypersensitive to mannitol, 
suggesting that the hypersensitivity of the transgenic lines to 
and Na^ ions are not due to general osmotic stress (Fig. 6D). 
The result, therefore, is consistent with the notion that overex- 
pression of the AtGluRl gene causes Ca^"^ deficiency. 

Cc^^ suppresses AtGluR2 -induced hypersensitivity to ionic 
stresses 

We, then, tested the effect of Ca^^ supplementation on the 
ionic stress response. Plants were grown in the salt stress media 
with supplementary Ca^^ for 2 weeks. As shown in Fig. 7, sup- 
plementing the media with 6mM CaClj dramatically attenu- 
ated the ion sensitivities of the y4/G/u/?2-overexpressing plants. 
It was observed that supplementation of Ca^"^ in the growth 
medium causes increase of the calcium contents in the plant 
shoots (Fig. 5A). Thus, supplementation of Ca^"^ in the medium 
in the presence of other ions should have increased relative 
concentration of Ca^"^ in the plant shoots. This, in turn, should 
have attenuated the effect of the cations in the transgenic lines. 
This interpretation is consistent with the notion that the trans- 
genic lines have lowered efficiency in Ca^^ utilization, mediat- 
ing the hypersensitivity to ionic stresses that can be rescued by 
increased Ca^"^ concentration. 

AtGluR2 expresses in vascular tissues 

In initial experiments, we tried to determine the distribu- 
tion of AtGluR2 mRNA among different organs by RNA gel 
blot analysis. However, we failed to detect any signal for 
AtGluR2 mRNA even when we used 25 ^g of total RNA (data 
not shown). However, RT-PCR analysis showed that the 
AtGluR2 gene is expressed in leaves, stems, flowers, and roots 
(data not shown). A more detailed localization of the AtGluR2 
gene expression was investigated using transgenic plants carry- 
ing the coding sequence of the E, coli p-glucuronidase (GUS) 
gene under the control of the AtGluR2 promoter. Activity of the 
GUS reporter gene was analyzed with the T2 progenies of three 
independent transgenic lines. The GUS activity was detected in 
the vascular tissues of both shoots and roots. Histochemical 
analysis in cross-sections revealed that the GUS activity was 
particularly strongly localized in cells adjacent to the conduct- 
ing vessels (Fig. 8). 
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Fig. 6 Sensitivity of transgenic plants to ionic stresses and osmotic stress. Plants were grown in l/4x B5 media supplemented with various con- 
centrations of KCl, NaCl, MgClj, or mannitol. The fresh weight is presented as a percentage relative to the fresh weight in I/4x B5 medium. Error 
bars represent the standard error of three independent experiments (25 seedlings were used for each measurement). (A) Sensitivity to KCl. (B) 
Sensitivity to NaCl. (C) Sensitivity to MgClj. (D) Sensitivity to mannitol. 



Discussion 

As an effort to understand the role of Ca^"*" in nutrition and 
the mechanisms underlying Ca^'*'-mediated adaptation to stress 
in plants, we have focused on the Arabidopsis homologs of 
ionotropic glutamate receptors. We characterized one member 
of this gene family as a candidate for a Ca^^ permeable chan- 
nel gene in higher plants. The family includes 19 genes that 
encode proteins with homology to manrunalian ionotropic 
glutamate receptors (Table 1). Ionotropic glutamate receptors 
(iGIuRs) are nonselective cation channels that are permeable to 
Ca^*^, K"*^, and Na"*" but not to Mg^^. However, the primary phys- 
iological function of these channels is thought to be regulation 
of Ca^"*" influx into cells. The Arabidopsis homologs show rela- 
tively low overall sequence identities (21% or less) with the 
mammalian iGluRs. However, the transmembrane domains 
(62% identity in TM III region) and the pore-forming region 
display a substantially high similarity; thus it is reasonable to 
presume that these plant proteins also function primarily as ion 
channels. The hydropathy plot of AtGluR2 supports the three- 
plus-one transmembrane topology as has been suggested in 
iCluRs. 

To investigate the physiological roles of the AtGluRl gene 
in plants, we generated transgenic lines of Arabidopsis that 



overexpress the AtGluR2 mRNA. The plants overexpressing 
AtGluR2 gene displayed symptoms that are characteristic indi- 
cators of Ca^^ deficiency in plant nutrition. The phenotype was 
exhibited in 19 independent transgenic lines and is therefore 
causally related to the introduced genes. The AtGluR2 gene 
overexpressing lines were also hypersensitive notably to 
and Na"^ ions. The ability to suppress the growth defects and 
ion sensitivity of the transgenic plants with exogenous Ca^"^ 
strongly suggests that these plants are deficient in Ca^^. 

Ca^^ deficiency can be caused by reduced-uptake of Ca^^ 
in the root, reduced transportation of Ca^"^ from root to shoot, or 
by reduced utilization of the transported Ca^^. Our data sug- 
gests that the major reason for Ca^"^ deficiency in the trans- 
genic lines is due to reduced utilization of the transported Ca^"*^. 

Developing a hypothesis that explains how overexpression 
of an ion channel protein might cause the reduced utilization of 
Ca^^ in transgenic plants is challenging. It is possible that the 
ectopic expression of the AtGluRl gene cause excess or Na"^ 
uptake if AtGluR2 is permeable to Ca^^ as well as monovalent 
cations as in the case of the mammalian iGluRs. Increased 
uptake of or Na*^ could compete with Ca^^, rendering 
reduced efficiency of Ca^^ utilization. An alternative possibil- 
ity would be that the overexpression of the AtGluRl gene 
might disturb the Ca^^ transport directly. If Arabidopsis gluta- 
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Fig. 7 The effects of supplemented Ca^* on ion hypersensitivity. (A) 
Plants growth on ionic stress conditions in the presence or absence of 
6 mM Ca^^ supplementation. Vector control plants are shown on left 
column and the transgenic A- 12 plants on the right column. The plants 
were grown for 2 weeks.(B) Relative fresh weight of the vector con- 
trol and the transgenic (A- 12, C-1) seedlings. Data are given as the 
mean±S£ of three independent experiments (25 seedlings were used 
for each measurement). Shown are values relative to the fresh weight 
of seedlings grown in l/4x B5 medium. 




Fig. 8 Expression pattern of the AtGluR2 gene in Arabidopsis. The 
GUS activity in transgenic Arabidopsis lines expressing the GUS 
reporter gene under the control of the AtGluR2 gene promoter was 
examined. The GUS activity was detected in the vascular tissues of the 
shoots (A) and the roots (C, D). Cross-section analysis of the GUS 
expression (B). The cross section was 7 ^m thick. P, phloem; X, 
xylem. Bar: 100 \xm. 



mate receptor homologs form heteromeric channels as in the 
case of the mammalian iOluRs, overexpression of a subunit 
could induce changes in overall channel activities. In the case 
that expression of the AtGluR2 protein results in heteromeric 
channels that are less efficient in Ca^^ transport, then overex- 
pression of AtGIuR2 might cause Ca^^ deficiency in low Ca^"^ 
nutrition. In mammalian cells, Ca^^ influx evoked by gluta- 
mate treatments is reduced by co-expression of die GluR2 sub- 
unit (Hollmann et al. 1991). 

The processes involved in unloading Ca^"^ from xylem 
vessels to the apoplast and in uptaking Ca^"^ by leaf cells are not 
fully understood. In fact, none of the mechanisms that regulate 
fluxes of Ca^" between various tissues of higher plants are well 
characterized. However, it is reasonable to assume that Ca^^ 
channels should be involved at certain steps to mediate Ca^* 
influx through plasma membranes. The promoter of the 
AtGluR2 gene is expressed along the vascular tissues, particu- 
.,20 larly in cells adjacent to conducting vessels, where Ca^^ is 



removed from the xylem sap. The characteristics of the 
deduced peptide sequence and the preferential expression of 
ih^AtGluRl gene in vascular tissues suggests that it may play a 
functional role in Ca^"^ unloading and allocation among differ- 
ent tissues. 

In any event, the AtGluRl transgenic plants are of interest 
as a starting point for understanding the molecular components 
for regulating efficiency of Ca^"^ utilization and for further 
functional characterization of the family of glutamate receptor 
channels in plants. In future studies, it will be of considerable 
interest to employ various heterologous expression systems to 
gain further insight into the biochemical and electrophysiologi- 
cal functions of these proteins. 
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It is well accepted that endogenous and environ- 
mental signals can influence cellular activities by 
changing [Ca^'**]^^ (Malh6 et al., 1998; Sanders et aL, 
1999). Despite the importance of this mechanism for 
coupling stimuli to responses (Malh6 et aL, 1998), the 
molecules responsible for generating increases in 
[Ca^'*']cyt during cell signaling in plants are not 
known at the genetic level. The results presented 
here raise the possibility that ligand-gated ion chan- 
nels in plants such as those predicted by the dis- 
covery of ionotropic glutamate receptor (iGluR)- 
like genes in Arabidopsis (Lam et al., 1998) are key 
components of a Ca^"^ influx mechanism important to 
signal transduction. 

In animal brains iGluR channels mediate fast chem- 
ical transmission across synapses by increasing the 
permeability of the post-synaptic cell membrane to 
K"^, Na"^, and Ca^"^ after binding Glu released by the 
presynaptic cell (Hille, 1992; Hollmann and Heine- 
mann, 1994; Dingledine et al., 1999). The resulting 
Ca^"*" entry in particular has been associated with 
long-term potentiation of the synapse, a physiologi- 
cal adaptation important to the learning process 
(Baudry and Lynch, 1993; Bliss and CoUinridge, 
1993). The Glu receptor homologs recently identified 
in plants (Lam et al., 1998) are too divergent from 
animal iGluRs to know with any certainty what li- 
gand(s) gate them, what ions are conducted in the 
open state, and in which membrane(s) of the cell they 
function (Chiu et aL, 1999). Thus the identification of 
iGluR sequences in the Arabidopsis genome raises 
intriguing questions about the physiological func- 
tions of neurotransmitter-gated chaimels in plant 
cells. 

The possibility that Glu gates Ca^^-permeable 
channels at the plasma membrane of plant cells was 
explored by measuring [Ca^^^J^^j in transgenic seed- 
lings expressing aequorin, a Ca^^ -sensitive lumin- 
escent protein (Knight et al., 1991). As shown in 
Figure lA, Glu application immediately triggered a 
very large, transient spike in [Ca^"^]cyt- ^ separate 
experiments the effect of Glu on membrane potential 
(^m) measured with intracellular microelec- 

trodes inserted into root apices of intact seedlings. 
Figure lA also shows that switching the bathing 
medium from 1 mM KCl to 1 mM K-Glu induced a 
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large and rapid depolarization of the membrane, as 
would be expected if the abrupt increase in [Ca^'^l^^t 
was due to Glu opening Ca^"^-permeable channels at 
the plasma membrane. The average peak change in 
induced by 1 mM Glu was 55 ± 7 mV (n = 6 
seedlings). This positive shift in V^, though consis- 
tent with Glu gating an inward electrogenic Ca^^ 
current across Sie plasma membrane, may also be 
due to secondary effects of the increased [Ca^"*']^^ on 
other ion transporters. Another scenario to consider 
is that Glu directly gates channels permeable to ions 
such as Cr (Cully et al., 1994) in addition to Ca^"^- 
permeable chaimels to cause the depolarization. And 
last, an electrogenic Glu-uptake mechanism (Boorer 
et al., 1996) may also contribute to the electrical re- 
sponse. Because these and perhaps other scenarios are 
not mutually exclusive, more electrophysiological 
studies of the connection between the large Glu-gated 
changes in [Ca^^l^y^ and the effect on are 
warranted. 

Figure IB shows that the magnitude and time 
course of the rapid increase in [Ca^'^]cyt was similar 
to the well-studied response to cold shock, i.e. in 
the micromolar concentration range and completed 
within several seconds (Knight et al., 1991, 1996; 
Lewis et aL, 1997). Treatment with 1 mM Glu induced 
a response that was tj^ically himdreds of fold higher 
than the control injection of equimolar KCl, which 
produced a touch response that may reflect Ca^^ 
entering the cytoplasm from internal stores (Haley et 
al., 1995; Legu6 et aL, 1997). The post-peak shoulder 
apparent in the selected response to Glu was often, 
but not always observed. Activation of iGluRs by Glu 
causes very similar Ca^"*^ changes in cells of the mam- 
malian nervous system (Kirischuk et aL, 1999; Obri- 
etan and van den Pol, 1999). 

If the increase in [Ca^'^'J^y^ triggered by Glu resulted 
at least in part from flux across the plasma membrane 
from the apoplasm, impermeant channel blockers 
and external chelators of Ca^"^ should reduce the 
response. The results in Figure 2A demonstrate that 
pretreatment with La'^'^, a frequently used blocker of 
plasma membrane Ca^^ channels, inhibited the Ca^"*" 
spike to the low level induced by the control treat- 
ment. Chelating extracellular Ca^"^ by pre-treating 
seedlings with EGTA was similarly inhibitory (Fig. 
2B). The combined evidence support our suggestion 
that Glu triggers an influx of Ca^^ across the plasma 
membrane and this leads to a dramatic change in 
[Ca^^lcyf Calcium-induced Ca^"^ -release from inter- 
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cific Ca^'^-charmel blocker (Lewis and Spalding, 
1998) and it may prevent the depolarization by block- 
ing a separate Glu-gated conductance in addition to 
the Ca^"^ pathway. Thus despite the fact that La^"^ 
blocks the Ca^^ flux and the membrane depolariza- 
tion, the exact relationship between the two Glu- 
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Figure 1. Clu triggers a large transient increase in [Ca^'^l^yj and an 
accompanying membrane depolarization. A, The red trace shows 
Ca^"^ -dependent luminescence from whole aequorin-expressing Ara- 
bidopsis seedlings (5- to 8-d-old) measured with a luminometer as 
described previously (Lewis et al., 1997), The black trace shows the 
response of V^, measured by impaling a cell near the root apex with 
an intracellular microelectrode as previously described (Spalding et 
al., 1999). Intact seedlings between 7- and 14-d-oId were used for 
the measurements. GIu at a final concentration of 1 mM was 
delivered as the salt. The pH was buffered at 5.7 with 2.3 mM MES 
[2-(N-morphoIino)-ethanesuIfonic acid]. 8, The kinetics and magni- 
tude of the change in ICa^'^Jcyt induced by Glu resembles the re- 
sponse to cold shock, but is much larger than the touch response 
induced by the control treatment. Cold shock was achieved by 
injecting 0°C 1 mM KCI into the luminometer cuvette, whereas the 
control treatment was room temperature 1 mM KCI. Glu was deliv- 
ered as 1 mM K-GIu and all solutions were buffered at pH 5.7. Arrow 
indicates the time of treatment. 



nal stores such as the vacuole may also contribute 
(Allen et al., 1995). 

La^"*" also blocked the depolarization triggered by 
Glu without affecting the resting V^^^ (Fig. 2C), indi- 
cating that the depolarization is a consequence of the 
inward Ca^*^ movement. However, La^"*" is not a spe- 
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Figure 2. Inhibitory effects of La^*^ and ECTA. A, Pretreatment of 
seedlings with 5 mM La^'^, an extracellular Ca^ "^-channel blocker, 
inhibited the change in ICa^^J^yt induced by Glu, but had a much 
lesser effect on the 100-fold smaller response to the control treat- 
ment. B, Chelating extracellular Ca^"^ by pretreatment with EGTA 
inhibited the Glu-induced Ca^"*" response. The data plotted are 
means (±sem) from three or four independent trials. C, Treatment of 
roots with La-*"*" blocked the Glu-induced depolarization without 
affecting the resting V^. Roots were treated with 5 mM LaClj before 
and during exposure to 1 mM K-Glu. 
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gated phenomena should be considered an open 
question. An alternative test would be to determine 
if EGTA treatment also inhibits the depolarization 
triggered by Glu, but stable recordings of are 
difficult to obtain when extracellular Ca^^ is de- 
pleted to an extent that significantly affects its avail- 
ability for inward fluxes. Patch-clamp studies of the 
ionic currents activated by Glu would be the pre- 
ferred means of obtaining a biophysical description 
of the depolarization mechanism. 

Glu is the primary natural ligand of iGluRs in the 
central nervous system although other non-native 
ligands are effective and have been used to classify 
receptor subtypes. The effectiveness of different li- 
gands was tested using the aequorin reporter plants. 
Figure 3 demonstrates that Glu was clearly the most 
effective agonist tested (note the logarithmic y-scale). 
a-amino-3-hydroxy-5-methylisoxazole-4-propionate 
and N-methyl-D-aspartate, potent agonists of animal 
iGluRs, did not induce a response above the control 
treatment (approximately 2% of the l-G1u response). 
These non-native agonists of animal iGluRs also did 
not activate the Synechocystis GLUO (Chen et al., 1999). 
It may be that affinity for a-amino-S-hydroxy-S- 
methylisoxazole-4-propionate and N-methyl-D-aspar- 
tate evolved in Glu receptors after the divergence of 
plants and animals. An alternative possibility is that 
Glu-gated Ca^"^ entry in Arabidopsis does not involve 
iGluR-like molecules, but instead some imrelated 
Ca^'^-permeable pathway lacking affinity for typical 
iGluR agonists is responsible for the phenomenon. 




Figure 3. Relative effectiveness of related compounds. l-CIu was 
much more effective than other potential agonists, including d-CIu 
and the animal iCIuR agonists, NMDA and AMPA. Note the loga- 
rithmic scale of the y axis, and that the response magnitudes are 
shown relative to the response induced by l-CIu. All compounds 
were administered at a final concentration of 1 mM. The plotted 
values are the means (±SE) of six independent trials. 
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The fact that the response to d-G1u was less than 
10% of the L-Glu response indicates high stereo- 
chemical specificity of the binding site(s) on what- 
ever molecules are responsible. Although Glu is 
clearly an effective ligand in a plausible concentra- 
tion range, other ligands may be more physiologi- 
cally important. The Arabidopsis genome contains 
several iGluR-like genes (Chiu et al., 1999) and that 
diversity may be matched by a similar diversity of 
agonists. 

Information on the effective concentration range of 
Glu would help to establish a physiological context 
for this ligand-gated response in plants. The change 
in [Ca^^](yt induced by Glu increased between 0.3 
and 3 mM, with the concentration for half-maximal 
response (EC50) being approximately 1 mM (data not 
shown). This value is approximately 10-fold greater 
than the typical value for prokaryotic and animal 
iGluRs, but very similar to the EC50 of CI "-permeable 
iGluRs from nematodes (Cully et al., 1994). 

If Glu or some other related small organic acid is 
the primary endogenous ligand, then it is important 
to consider how and when the external ligand- 
binding site would experience 0.3 to 3 mM concen- 
trations. Anion channels at the plasma membrane of 
plant cells are known to function in the transduction 
of several signals important to plant growth and 
development (Ward et al., 1995). These channels are 
relatively non-selective among anions and may con- 
duct significant efflux of dicarboxylic anions such as 
malate (Hedrich, 1994; Schmidt and Schroeder, 1994), 
and therefore perhaps Glu, as well. When environ- 
mental or endogenous signals activate such aruon 
channels, apoplastic Glu concentration may rise into 
the effective range, causing a transient change in 
[Ca^"^]cyt that serves to couple a stimulus to down- 
stream responses. This hypothetical scenario may be 
most plausible in roots, where anion-channel medi- 
ated release of dicarboxylic acids has been proposed 
as a mechanism for combating Al^^ stress (Delhaize 
and Ryan, 1995; Ryan et al., 1997). Perhaps it is no 
coincidence that dissection experiments revealed 
most of the Ca^"*^ signal recorded from intact Arabi- 
dopsis seedlings was contributed by the root; leaves 
and cotyledons of yotmg plants displayed smaller 
Glu responses (data not shown). 

The results presented here form the basis of our 
proposition that a key element of a mechanism for 
altering [Ca^'^^lcyt in plant cells during signaling is 
similar to that responsible for neurotransmitter ac- 
tion in the central nervous system of animals. The 
evidence would be bolstered considerably if muta- 
tional studies revealed a link between specific iGluR- 
like genes and Glu-gated Ca^"^ fluxes. Plant biologists 
interested in Ca^^ signaling are presently particu- 
larly well equipped to test this connection because 
there is a wealth of published details about iGluR- 
mediated Ca^^ signaling in neurons, the Arabidopsis 
genome is essentially sequenced and searchable, so- 
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phisticated reverse-genetic strategies are very practi- 
cal, and electrophysiological techniques can measure 
function with high resolution. The stage for exciting 
developments in Ca^"^ signaling is set. 
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